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Abstract 

The astronomical dark matter could be made of weakly interacting and massive 
particles. If so, these species would be abundant inside the Milky Way, where 
they would continuously annihilate and produce cosmic rays. Those annihilation 
products are potentially detectable at the Earth, and could provide indirect 
clues for the presence of dark matter species within the Galaxy. We will review 
here the various cosmic radiations which the dark matter can produce. We will 
examine how they propagate throughout the Milky Way and compare the dark 
matter yields with what pure astrophysical processes are expected to generate. 
The presence of dark matter substructures might enhance the signals and will 
be briefly discussed. 

Resume 

La matiere noire astronomique pourrait etre constitute de particules mas- 
sives aux interactions evanescentes. Si tel etait le cas, ces particules se re- 
trouveraient en abondance au sein de la Voie Lactee ou elles s'annihileraient 
en permanence, produisant de multiples radiations cosmiques. Celles-ci sont 
eventucllcmcnt visibles de la Terre et constituent des lors des sortes d'empreintes 
spectrales, veritables signatures indirectes des candidats potentiels a la matiere 
noire galactique. Dans cet article, nous passons en revue les differentes especes 
cosmiques susceptibles d'etre produites, et comparons leur flux avec celui des 
radiations engendrees par les processus astrophysiques conventionnels. L'exis- 
tence de condensations de matiere noire est brievement discutee. Le taux d'an- 
nihilation pourrait etre amplifie au sein de telles structures, conduisant a des 
signatures indirectes plus intenses que dans le cas d'un halo galactique lisse. 
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1. The messengers of dark matter annihilation 

The nature of the astronomical dark matter (DM) is still unknown. This 
component, which contributes a quarter to the energy balance of the universe, 
cannot be made of atoms and nuclei like ordinary matter. An exciting, and quite 
plausible possibility, lies in the existence of a population of weakly interacting 
and massive particles, dubbed the acronym WIMP. Should DM species pervade 
the halo of the Milky Way, their mutual annihilations would yield several indirect 
signatures. These are potentially detectable at the Earth under the form of 
spectral distortions appearing in various cosmic radiations 

X + X^ qq,W + W-,...^p,D : e+ 1 kv's . (1) 

Detection of the DM annihilation products has motivated the spectacular devel- 
opment of several new experimental techniques. Searches for antiprotons and 
positrons are performed by balloon and satellite-borne devices. Because the 
flux depends on the square of the WIMP density n x , the limit which may be 
set on the annihilation cross section scales very roughly as That type of 
search is a priori mostly sensitive to small WIMP masses. Actually, the recent 
discovery of a positron excess above 10 GeV by the PAMELA collaboration has 
triggered a lot of excitement in the field, and many new possibilities have been 
explored, which may lead to a strong signal even in the TeV window. Alas, the 
DM interpretation of the positron excess raises more problems than it solves, 
and the current explanation of that anomaly lies in the existence of nearby 
pulsars. High-energy photons are detected both by air Chercnkov telescopes 
(ACT) and by satellite-borne instruments. The WIMP annihilation rate and 
hence the gamma- ray signal both scale as m~ 2 . Because of the background in 
which that signal is swamped, the experimental reach on the annihilation cross 
section approximately scales as the mass m x . A smoking gun signature of the 
presence of DM particles would be the detection of a hot spot in the gamma-ray 
sky. The neutrino channel is mostly sensitive to large values of m x . The limit 
which may be set on the annihilation cross section does not depend too much 
on the WIMP mass and this channel is complementary to the other searches. 
Large areas under the ice cap of the South Pole are equipped for the detection of 
high-energy up-going muons, and the theoretical predictions start to be within 
reach and will soon be checked. 

Various species can be produced by WIMP annihilations among which are 
antimatter cosmic rays, high-energy photons and neutrinos. The corresponding 
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rate <?dm (x, E) for the production of these particles depends on their energy E 
and is related to the WIMP annihilation cross section cr ann through 

g DM (x,£) = V (cw) ( — } f(E) . (2) 

The coefficient rj is a quantum factor equal to 1/2 for a self-conjugate particle 
like a Majorana fermion or to 1/4 otherwise. The annihilation cross section 
is averaged over the momenta of the incoming DM particles to yield (cr ann w), 
whose value depends on the specific microscopic interactions involved in the an- 
nihilation process. The DM density at location x is denoted by p(x), while f(E) 
stands generically for the energy distribution dN/dE of the species generated 
in a single annihilation event. 

Once produced, the cosmic radiations propagate within the Milky Way and 
eventually reach the Earth. A key issue is the transport of charged cosmic rays 
throughout the magnetic fields of the Galaxy. An overview of the subject is given 
in section[2j This allows us to focus on antimatter particles. Although produced 
by conventional astrophysical processes, these species are not abundant. They 
are promising targets insofar as their spectra could be distorted should addi- 
tional (and exotic) sources, of DM origin in our case, operate. Sections |3j [4] and 
[5] are respectively devoted to antiprotons, positrons and antideuterons. Atten- 
tion is paid to the astrophysical backgrounds inside which the DM signals are 
buried. The PAMELA positron excess is also discussed, and its interpretation 
in terms of WIMPs is shown to have difficulties. The gamma-ray and neutrino 
skies are the subjects of sections [6] andj7l Finally, the possibility that DM is 
clumpy is rapidly examined in section |4.3| The existence of substructures has 
been extensively used in the literature as a pretext for arbitrarily enhancing 
the fluxes from DM origin. These are generally depressingly weaker than the 
backgrounds inside which they are hidden. However, a recent analysis has es- 
tablished that the boost factor from DM dumpiness cannot exceed on average 
a generous factor of ~ 20. Alternatively nothing precludes a DM clump to lie 
close to the Earth and produce locally an intense flux of cosmic radiations. But 
the odds for this to happen are not large. 

For a more extensive view of the dark matter enigma, we refer the reader to 
the complementary contributions of Silk [1 ] , Hooper & Tait [5] , and Armengaud 
[3] in the present issue. 

2. An overview of cosmic ray transport 

We will focus this section on the transport of charged cosmic rays. Photons 
and neutrinos propagate along straight lines, and their fluxes at a distance r from 
a point source decrease as r -2 . Some extinction can take place along the line of 
sight though, but this is negligible for gamma-rays in the GeV to TeV energy 
range. The situation for neutrinos is more complicated as oscillations among 
the various families is known to occur in empty space, and is even resonantly 
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enhanced inside the Sun and the Earth. Neutrino transport is broached in 
section 17.11 

Once charged cosmic rays have been accelerated, or created by spallations 
of other nuclear species or through DM annihilation, they propagate through 
the Galactic magnetic field and are deflected by its irregularities : the Alfven 
waves. In the regime where the magnetic turbulence is strong, which is actually 
the case for the Milky Way, cosmic ray (CR) transport needs to be investigated 
numerically. Monte Carlo simulations [4] indicate that it is similar to space 
diffusion with a coefficient 

K(E) = K Q (3 (K/l GV) S , (3) 

which increases as a power law with the rigidity 1Z = p/q of the particle. In 
addition, because the scattering centers drift inside the Milky Way with a veloc- 
ity V a ~ 20 to 100 km s _1 , a second order Fermi mechanism is responsible for 
some mild diffusive reacceleration. Its coefficient Kee depends on the particle 
velocity /? and total energy E, and is related to the space diffusion coefficient 
K{E) through 

K EE = -V 2 (A) 
EE 9 a K{E) ■ { ' 

In the case of positrons, diffusive reacceleration is completely dominated by 
energy losses. Finally, Galactic convection wipes cosmic rays away from the 
disk with a velocity Vc ~ 5 to 15 km s _1 . 

Radio observations indicate that the magnetic fields of galaxies are not con- 
fined to the visible matter which they host, but extend far away in space. An 
illustration is provided by the maps of the radio continuum halo of NGC 4631, 
taken at 610 and 1412 MHz The synchrotron emission from electrons spiral- 
ing inside the magnetic field of NGC 4631 probes the structure of its magnetic 
halo. This halo, inside which charged particles are trapped and diffuse, extends 
well above and beneath the luminous disk. Inspired by this example, we can 
model the magnetic halo of our Galaxy (the so-called diffusive halo or DH) as 
a thick disk which matches the circular structure of the Milk Way, as shown 
in Fig. [I] The Galactic disk of stars and gas, where primary cosmic rays are 
accelerated, lies in the middle. Primary species, such as CR protons, helium 
nuclei and electrons, are believed to be accelerated by the shock waves driven 
by supernova explosions. These take place mostly in the Galactic disk which 
extends radially 20 kpc from its center, and has a half-thickness h of 100 pc. 
Confinement layers, where cosmic rays are trapped by diffusion, lie above and 
beneath this thin disk of gas. The intergalactic medium starts at the vertical 
boundaries z — ±L, as well as beyond a radius of r = R = 20 kpc. Notice that 
the half-thickness L of the diffusive halo is not known and reasonable values 
range from 1 to 15 kpc. The diffusion coefficient K is assumed to be the same 
everywhere inside the DH, whereas the convective velocity is exclusively vertical 
with component Vc(z) = Vc sign(z). This Galactic wind, which is produced by 
the bulk of the disk stars like the Sun, drifts away from its progenitors along 
the vertical directions, hence the particular form assumed here for Vc- 
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Figure 1: Schematic edge-on view of the Milky Way diffusive halo (DH) as seen by a cosmic 
ray physicist. The stellar and gaseous disk is sandwiched between two thick layers which 
contain turbulent magnetic fields. After having been accelerated by supernova driven shock 
waves or produced by DM species annihilating in the Galactic halo, cosmic rays diffuse on 
magnetic inhomogeneities and are wiped away by a Galactic wind with velocity Vc- They can 
lose energy and are also mildly subject to diffusive reacceleration. The former process is by 
far the dominant one in the case of electrons and positrons. This diagram has been borrowed 
from the review [5]. 



The master equation for CR transport may be expressed as 
^ + d z (V c iP) - KA^ + d E {b loss (E)ij - K EE {E)d E ^} = q( X ,E) , (5) 

where ip = dn/dE denotes the CR space and energy distribution function. This 
relation applies to any species (protons, antiprotons or positrons) as long as 
the rates for production q and energy loss b loss (E) are properly accounted for. 
Most of the analyses devoted to Galactic cosmic rays are based on the assump- 
tion that the acceleration and propagation of charged particles have reached 
a steady state. The diffusion time of a 10 GeV proton throughout the mag- 
netic inhomogeneities of the Milky Way disk is of order 7 million years whereas 
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the rate of supernova explosions is 1 to 3 every century. The case for steady 
state is made even stronger when the escape time ~ L 2 /K from the DH as a 
whole is considered instead of the Galactic disk residence time. A 10 GeV proton 
spends on average 2.9 x 10 s years before it escapes into intergalactic space. Dur- 
ing that period, approximately 3 to 9 million supernova explosions have taken 
place. It is then reasonable to describe these sources as if they were continuously 
spread along the Galactic disk and were steadily accelerating charged particles. 
Different methods have been proposed to solve the CR diffusion equation j|5]). 
The purely numerical approach followed in the GALPROP public code 7jFTis 
based on the well-known Crank-Nicholson semi-implicit scheme, where the time- 
dependent equation is evolved on a space and energy grid until convergence is 
reached. This method allows to treat inhomogeneous cases where, for instance, 
K or Vc depend on the location x within the DH. It is time consuming though, 
and particular attention needs to be paid to the boundary conditions. More 
recently, the DRAGON package [8 ^ has been made publicly available. Most of 
the building blocks of DRAGON come from GALPROP, but this code relies on 
a faster solver. 

Another route has been taken by the USINE collaboration, based on two 
semi-analytic methods. 

(i) The Bessel expansion technique takes advantage of the axial symmetry of the 
DH and enforces a vanishing cosmic ray flux at a distance R = 20 kpc from the 
rotation axis of the Galaxy. This condition is actually implemented naturally 
by the following series expansion for ip 

+oo 

1>(x,E) = i/>(r,z,E) = Y / ^(^ E )M^r/R) . (6) 
i=l 

The Bessel function of zeroth order J vanishes at the points Oj. The radial 
dependence of tp is now taken into account by the set of its Bessel transforms 
Pi(z, E). The source term q(x, E) = q(r, z, E) may also be Bessel expanded into 
the corresponding functions Qi(z,E) so that the master equation ^ becomes 

d z {VcP l )-KdlP l + K\^ 2 p i + d E {b loss {E)P l -K EE {E)d E P l } = Q t (z,E) 

(7) 

In the case of CR nuclei and antinuclei, energy losses and diffusive reacceleration 
are confined inside the Galactic disk, considered here as infinitely thin, and an 
effective term 2 h S(z) should be in factor of the energy derivative. 

(ii) The solution of the master equation ^ may also be generically expressed 
as the integral 

^(x ,£) = fdEsf d 3 x s G( XQ ,E^x s ,E s )q( Xs ,E s ) . (8) 

J J DH 



1. See the web site 

2. See the web site 



http://galprop.stanford.edu/ and references therein 



http://www.desy.de/~maccione/DRAGON/ and references therein 
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The energy E$ at the source runs over a range which depends on the CR species 
and on the production mechanism. The space integral is performed over the 
DH. The convolution ^ involves the Green function G which describes the 
probability for a cosmic ray produced at location xg with energy Eg to reach, 
for instance, the Earth located at x , where it is detected with energy E. The 
cosmic ray space and energy density tp can be translated into the differential flux 
$ = (j3/ An) ip where f3 stands for the particle velocity. This flux is expressed in 
units of m -2 s~ x sr _1 GeV -1 . The flux of antinuclei or positrons produced by 
WIMP annihilations may be written as the product 

*(©, E ) =jr f dE s f(E s ) I(E, Eg) , (9) 

where the information related to particle physics has been factored out in 

T = V~KnnV) ( — V • (10) 
4tt [m x ) 

The energy distribution f(Es) describes the spectrum at the source and depends 
on the details of the WIMP annihilation mechanism. The information on the 
Galactic DM density profile p, as well as on the propagation of cosmic rays 
within the Milky Way DH, is summarized in the halo integral 

I(E,E S ) = ( d 3 x s G(x e ,E^x s ,E s ) (^H , (ii) 

where the solar neighborhood DM density is denoted by p & . The halo integral 
I(E, Es) is a key ingredient for the derivation of the flux at the Earth of an 
antimatter species produced inside the Galactic DH by WIMP annihilations. 
The spatial reach of the Green function G depends on the nature of the CR 
particle (antinuclei or positrons) and on the energies E and Es- This range 
delineates the region of the Milky Way from which most of the signal detected 
at the Earth originates. It corresponds to the extension of the so-called horizon 
beyond which the Green function vanishes. The antiproton horizon reaches the 
Galactic center at high energy, and the antiproton flux at the Earth starts to be 
sensitive to the DM distribution there. On the contrary, the positron horizon 
shrinks at high energy. Above tens of GeV, positrons detected at the Earth 
originate from its vicinity. 

The normalization coefficient Kq, the spectral index 6, the Galactic drift 
velocity Vc and the Alfven velocity V a are all unknown. We have already men- 
tioned that such is the case for the DH half-thickness L. This situation can be 
remedied with the help of the boron to carbon ratio B/C which is quite sensitive 
to CR transport and which may be used as a constraint. The three propagation 
models featured in table [T] have been drawn from [TU]. The MED configuration 
provides the best fit to the B/C measurements whereas the MIN and MAX 
models lead respectively to the minimal and maximal allowed antiproton fluxes 
which can be produced by WIMP annihilations. The three sets of parameters 
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Case 


5 


K [kpcVMyr] 


L [kpc] 


V c [km/s] 


V a [km/s] 


MIN 


0.85 


0.0016 


1 


13.5 


22.4 


MED 


0.70 


0.0112 


4 


12 


52.9 


MAX 


0.46 


0.0765 


15 


5 


117.6 



Table 1: Typical combinations of diffusion parameters that are compatible with the B/C 
analysis 3 . As shown in |10) . these propagation models correspond respectively to minimal, 
medium and maximal primary antiproton fluxes. 



of table [T] belong to an ensemble of more than 1,600 other models which have 
been shown in [9j to be also compatible with the B/C ratio. That analysis has 
been recently improved by several groups. In particular, the USINE collabora- 
tion has constrained the CR propagation parameters with the help of a Markov 
chain Monte Carlo technique [llj which allows to explore rapidly the parameter 
space. The resulting best fit is somewhat different from what has been derived 
ten years ago, but one of the important messages is the difficulty to get a precise 
value for L which still lies between 1 and 16 kpc. 

3. Antiprotons as a robust probe for DM species 

3.1. Calculation of the antiproton flux at the Earth 

Antiprotons are produced during the collisions undergone by primary CR 
nuclei on the interstellar gas, i.e., within the Galactic disk. Because they are not 
directly injected in the interstellar medium (ISM) but are sourced by primary 
species, these astrophysical antiprotons are dubbed secondaries. The rate with 
which they are produced may be expressed as 

f +oc da 
q s p cc (r, E p ) = n H x (3 p %{r, E p ) x dE p x—(E p -+ E p ) , (12) 

Je« d£j p 

in the case of interactions between CR protons and hydrogen atoms. The various 
contributions from the spallations of interstellar H and He by CR protons and 
alpha particles need to be taken into account. Details on how the cross sections 
of these processes are parameterized are given in [12 and [13 . In addition to 
this conventional mechanism, antiprotons may be directly produced as primary 
CR species by DM annihilation. The corresponding source term qp Ilm (r, z, E) 
has already been discussed and is generically given by expression ([2]), where 
f(E) stands here for the antiproton spectrum dN p /dE p . Notice that WIMP 
annihilations take place all over the diffusive halo and are not restricted to the 
z = region. We therefore anticipate a different sensitivity of that component 
to the CR propagation parameters than in the case of secondary antiprotons. 

Once produced, antiprotons propagate inside the DH. They can collide elasti- 
cally on interstellar H and He atoms. However, they are preferentially scattered 
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forward, so that these interactions are innocuous and will be disregarded. An- 
tiprotons can also annihilate on interstellar H and He. This leads to a negative 
source term — Tp nn ipp, where the annihilation rate Tp nn is defined as 

= /3p n H + /3p n Hc . (13) 

The annihilation cross section a^™ can be borrowed from P3HH] and multiplied 
by a factor of 4 2 / 3 ~ 2.5, taking into account the higher geometric cross section, 
to get o-p He- The average hydrogen rtjj and helium riHo densities in the Galactic 
disk can be respectively averaged to 0.9 and 0.1 cm" 3 . Last but not least, a 
tertiary component arises from the inelastic and non-annihilating interactions 
which antiprotons undergo with the ISM. Antiprotons can actually collide on a 
nucleon at rest and transfer enough energy to excite it as a A resonance. This 
mechanism redistributes antiprotons toward lower energies and flattens their 
spectrum as shown in |16j . It yields the source term 



P g pX (4 -> E p ) n H % ^(r, 4) dE P 

VpH^px(Ep) «H Pp i>p{r,Ep) , (14) 

where the inelastic and non-annihilating differential cross section in this expres- 
sion can be approximated by dapB_^px/dEp ~ <XpH->j5x/7p- This parameteri- 
zation can be improved by using the Anderson prescription |17j as described in 
[T5J[T5]. In relation (14 1, the initial antiproton kinetic energy is denoted by Tp. 
In order to take into account elastic scatterings on helium, one simply has to 
replace the hydrogen density by nu + 4 2/,3 riHe- 

The antiproton background and DM signal are reliably computed with the 
Bessel expansion approach. This method encodes directly the presence of radial 
boundaries for the magnetic halo. The full expression for the master equation 
describing the Bessel transformed antiproton distribution function Pi(z, E) may 
be expressed as 

d z {V c Pi) - KdlPi + K^Xpi 
+ 2hS(z)d E {b loss {E)P l - K BB {E)d E Pi} = 

-2^)1^^ + Qf™{z,E) + 2hS(z) {q;< ■ Q';;;} . (15) 

The various source terms have been specified above. Annihilations of antipro- 
tons take place inside the Galactic disk. This is also true for the production 
of secondary and tertiary antiprotons. All those processes involve the ISM. We 
have already mentioned, as a side remark on equation ([7]), that energy losses 
and diffusive reacceleration are also confined inside the Galactic disk. The ap- 
proximation of infinitely thinness for the disk leads to the effective term 2 h 5(z) 
in the above expression. Integrating it through the infinitely thin disk, along the 
vertical axis z, leads finally to a diffusion equation in energy which the Bessel 
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transforms Pj(0, E) fulfill 



AiPi(0,E) + 2hd E {b loss {E)P l {0,E) - K E E(E)d E P t {0,E)} = 

2h {QfS + QfZ} + 2 / dzQl™(z,E) e 2K Fl {z) . (16) 
Jo 

The coefficients Ai which appear in the above expression are given by 

ME) = V c + 2hT^ n (E) + K{E)S iC oth^^ 

where Sf = (V c /K) 2 + (2 a t /R) 2 , (17) 
while the vertical functions ^(z) are defined as 

F t {z) = sinh[|(L-z)]/sinh[|i] . (18) 



For each Bessel order i, the integro-differential equation (16) is solved following 
the procedure explained in the appendix B of |20j . Evolving in time an initial 
burst with the help of a semi-implicit Crank-Nicholson scheme allows to check 
the convergence of that procedure. The tertiary component is computed by re- 
injecting several times the total antiproton yield ipp in the integral over energy 



of relation (14 1. Convergence is reached very rapidly. 



3.2. Antiproton background and signal 

The spallations of interstellar H and He by cosmic ray primaries, essentially 
protons and alpha particles, produce an irreducible background of secondary 
antiprotons inside which the signal from DM species is hidden. The precise 
determination of this background is crucial in order to disentangle a possible 
WIMP signature. The semi-analytic treatment of CR propagation which has 
been discussed in section [2j and which is based on the Bessel expansion ([6]), is 
a convenient framework to derive the theoretical uncertainties associated to the 
various parameters at stake, namely K , S, V a , Vc and the DH half-thickness 
L. The space of these propagation parameters has been extensively scanned 
in order to select the allowed regions where the predictions on B/C, a typi- 
cal CR secondary to primary ratio, match the observations. Several hundreds 
of different propagation models have survived that test. The propagation pa- 
rameters are thus only loosely constrained by the CR nuclei abundances so far 
observed. The yellow band presented in each of the figures of this section is 
actually the envelope of the secondary antiproton spectra computed with the 
set of those ~ 1,600 different propagation models found in [9] to pass the B/C 
test. This band comprises the theoretical uncertainty in the determination of 
the secondary antiproton flux. It is confined by the MIN and MAX configu- 
rations of table [T] As a first observation, notice how narrow the uncertainty 
strip is between ~ 10 and 100 GeV. In spite of this, the theoretical predictions 
are in remarkable agreement with the recent antiproton measurements [24] of 
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Figure 2: The astrophysical antiproton background lies within the yellow band whose thick- 
ness indicates the theoretical uncertainty arising from CR propagation. This band actually 
encompasses the secondary antiproton yields computed for more than 1,600 different sets of 
CR propagation parameters, all compatible with the B/C ratio. The primary antiproton 
fluxes are produced by the annihilation of a Kaluza-Klein WIMP that appears in the higher 
dimensional warped Grand Unified Theories of j2 It 122} . The mass of that DM species (dubbed 
LZP) has been varied from 30 to 70 GeV with a Kaluza-Klein scale A/kk of 3 TeV. When 
the LZP mass is close to M z o/2, the annihilation becomes resonant and the primary signal 
may even exceed the background. The MAX diffusion parameters of table ^ have been used 
for the LZP antiproton spectra, with a canonical isothermal DM distribution. This figure has 
been borrowed from | 23| . 

the PAMELA collaboration. These span 4 orders of magnitude between 1 and 
100 GeV. At higher energies, the yellow band widens as a result of the energy 
dependence of the diffusion coefficient K. We expect antiproton propagation to 
be dominated by pure diffusion in that energy range. A very crude approxima- 
tion for the antiproton Green function is obtained in the limit where the DH is 
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infinite and may be expressed as 



where r ffi denotes the distance between the Earth and the source. The antipro- 
ton flux is expected to scale as $ p /K(E) oc E~ a ~ s , where a is the spectral 
index of the proton flux at high energies. From the B/C analysis, the spectral 
index 5 of the diffusion coefficient K may take any value between 0.46 and 0.85. 
Its spread = 0.4 translates for the antiproton flux into a factor of 3 of uncer- 
tainty at 1 TeV. Notice that PAMELA and AMS-02 will considerably improve 
the measurements of the CR nuclei abundances, with a determination of the 
B/C ratio to a better accuracy and over a wider energy range than available 
so far. This will translate into improved constraints on the propagation param- 
eters and eventually into a thinner uncertainty yellow strip for the antiproton 
astrophysical background. 

The antiproton signal from annihilating DM species leads to a primary com- 
ponent directly produced throughout the DH. It depends on many unknown 
ingredients as is clear from relation ([2]). The annihilation cross section (cr ann w) 
at freeze-out is related to the WIMP cosmological relic abundance through 

„ 3 x 10~ 27 cm 3 s- 1 

tt^h 2 ~ ; : . 20 

A value of f2 dm ~ 0.21 translates into a typical WIMP annihilation cross section 
of order 3 x 10 -26 cm 3 s . Although this is strictly true only at decoupling, 
we can keep in mind that value as a benchmark (see below though). The an- 
tiproton spectrum f(E) at the source is also model dependent. Antiprotons are 
generated through quark or gauge boson jets. Since the WIMPs are at rest in 
the Galactic frame, the antiproton spectrum tends to flatten at low energy. The 
production of tertiary antiprotons leading to the same effect, we do expect a flat 
spectrum below a few GeV, whatever the production mechanism. The primary 
and secondary fluxes displayed in Fig. [2] and [3] have actually fairly similar spec- 
tra at low energy. The WIMP mass is a key ingredient. The production rate 
9dm scales as m~ 2 . The largest DM antiproton signals are generally expected 
in the case of light DM candidates. For illustration, we have selected among 
the realm of possible DM candidates a particular Kaluza-Klein WIMP which 
appears in the context of higher dimensional warped Grand Unified Theories 
[2"T1 In these models, a stable KK fermion can arise as a consequence of 
imposing proton stability in a way very reminiscent of R-parity stabilizing the 
lightest supersymmetric particle in supersymmetric models. The symmetry is 
called Z3 and the Lightest Z3 Particle (LZP) is stable since it cannot decay into 
standard model particles. It is actually identified with a KK Dirac right-handed 
(RH) neutrino with a mass in the 1 GeV to 1 TeV range. This RH neutrino 
has gauge interactions in particular with additional KK Z' gauge bosons. Nev- 
ertheless, its interactions with ordinary matter are feeble because they involve 
heavy gauge bosons with a mass Mkk > 3 TeV. The primary antiproton flux at 



12 




° N 
o J 



-A3D 



VOJ.$ 







CO 
03 




o 


+ 


CO 


03 




o 


C3 


03 




CO 


N 


03 


CD 




03 








o 


03 


3 


u: 


CO 


5 






~i) 


CO 


Oh 


in 


s 


W 


PJ 


< 


S 


< 






O 


< 


a. 




■ 


• 


•4 


• 



> 

o 
o 

m 




> e 



6=r 



[l-A30 I- JS I- S Y0l# 



fa ft . 



13 



the top of the atmosphere (TOA) is plotted in Fig. [2] as a function of antiproton 
kinetic energy for five different values of the LZP mass. The most optimistic 
Galactic diffusion scheme MAX as well as a canonical isothermal Galactic DM 
halo have been assumed for the primary signal. The curves corresponding to 
Mlzp = 40 (short dashed magenta) and 50 GeV (long dashed red) exceed the 
background and should have already led to a detection would our assumptions 
on Galactic diffusion and halo profile be correct. For Mlzp = Mzo/2, the LZP 
annihilation is actually driven by the Z-resonance and is significantly enhanced. 

The antiproton DM signal is more sensitive to Galactic CR propagation than 
the secondary component. The latter is generated through the interactions of 
CR protons and helium nuclei with the ISM, a process quite similar to the pro- 
duction of boron through the spallation of carbon or nitrogen nuclei impinging 
on interstellar gas. Although the subnuclcar processes at stake are different, 
secondary antiprotons share a few similarities with boron nuclei. In particular, 
the production site is the Galactic disk in both cases. It is no surprise then if 
the yellow band of uncertainty is so narrow once the B/C constraint has been 
taken into account. Such is not the case for the antiproton DM signal which is 
produced all over the magnetic halo, in regions far away from the disk. We do 
expect then a larger variance for that signal than for the background as the CR 
propagation parameters are varied over the range allowed by the B/C measure- 
ments. In the left panel of Fig. [3j the primary antiproton flux produced by a 40 
or 50 GeV LZP species decreases by two orders of magnitude between the most 
optimistic (MAX) and the most pessimistic (MIN) diffusion cases of table [IJ In 
the last configuration, the antiproton signal is now well below the background. 

The distribution of DM inside galaxies is still an open question. From one 
side, results from cosmological N-body simulations in A-CDM models [251 023 
indicate a universal and coreless DM density profile. At small radii, the DM 
density diverges with the distance r from the Galactic center as r -7 , with 7 ~ 1 
to 1.5. Other results, also obtained from simulations of halo formation, strongly 
disfavour a singularity as steep as 1.5 and seem to point toward slopes logarith- 
mically dependent on the distance from the Galactic center and no steeper than 
~ 1.2 J27J HH]- On the other side, several analyses of rotational curves observed 
for galaxies of different morphological types [2H1 ISO] put serious doubts on the 
existence of DM cusps in the central regions of the considered objects. Instead 
of a central singularity, these studies rather suggest a cored DM distribution, 
flattened toward the central regions. The DM profile inside the Milky Way can 
be parameterized by the generic distribution 



where r = 8.5 kpc is the distance of the solar system from the Galactic center. 
The local DM density has been set equal to p© = 0.3 GeV cm -3 — note 
that this canonical value is slightly lower than more recent estimates made in 
e.g. Refs. [211 [321 133 GS] , which agree on a central value around 0.4 GeV/cm~ 3 . 
In the case of the pseudo-isothermal profile [35], the typical length scale a is 




(21) 
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the radius of the central core. The profile indices a, /3 and 7 for the DM 
distributions considered in Fig. [2] and [3] are indicated in table [2] As is clear in 



Halo model 


a 


P 


7 


a [kpc] 


Cored isothermal [35 


2 


2 





4 


Navarro, Frenk & White [25 


1 


3 


1 


25 


Moore [26] 


1.5 


3 


1.3 


30 



Table 2: Dark matter distribution profiles in the Milky Way. 



the right panel of Fig. [3j where the case of a 60 GeV LZP is considered, the 
Galactic DM distribution p(x) is also a source of uncertainty in the calculation of 
the antiproton DM signal. The MAX set of CR propagation parameters, which 
has been assumed in that example, makes it possible for primary antiprotons 
produced at the central Galactic cusp to reach the solar circle. The degeneracy 
among the various DM distributions is not lifted if the MIN propagation model 
replaces the MAX configuration. In that case, the three colored curves are one 
and the same. 

Because the astrophysical background of secondary antiprotons is well under 
control once the B/C constraint is taken into account, antiprotons offer a unique 
probe of the presence of annihilating DM species within the Galactic halo. This 
is especially true for light candidates. The DM signal is expected to be the 
strongest in that case. As a matter of fact, recent observations collected by the 
direct detection experiments CDMS-II [3(5] and CoGENT [37J are compatible 
with a light WIMP with mass around 10 GeV. Should such a particle couple 
preferentially to quarks, it would generate a sizable amount of antiprotons at 
low energy, and would be excluded by observations as shown in [38]. A light 
neutralino arises in supersymmetric extensions of the standard model with an 
extra gauge singlet superfield. Generic models embedding an additional Majo- 
rana fermion associated to two new scalar fields are constrained by antiprotons, 
although some regions of the parameter space still pass the test [35] ■ 

The annihilation rate of a heavy WIMP is suppressed by a factor of m~ 2 . 
The antiproton DM signal is hidden in the background, unless the annihilation 
cross section (er arm v) is abnormally large. This situation may happen in su- 
persymmetric models where the neutralino is almost a pure Wino, as expected 
for example in anomaly mediated supersymmetry breaking (AMSB) scenarios 
[40j . For Winos, the preferred mass from relic density requirements is peaked 
at about 1.7 TeV [41]. Non-perturbative binding energy effects then result in 
greatly enhanced annihilation cross sections today, when the neutralinos have 
very small Galactic velocities [HI 23J. In this limit, heavy Winos annihilate 
almost exclusively into gauge bosons. The case of a heavy Wino particle is 
featured in Fig. E] The annihilation cross section today (er ann i>) is equal to 
1.02 x 10 - 24 cm 3 !" 1 [13], well above the decoupling value. That species anni- 
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Figure 4: The yellow band features the expected antiproton background for the full range of 
diffusion parameters allowed by the B/C ratio. A heavy WIMP is also considered. This DM 
species is almost a pure Wino and its annihilation cross section is significantly enhanced today 
by non-perturbative, binding energy effects. The corresponding primary (long dashed) and 
total (solid) fluxes have been derived for a NFW halo profile and for the MED set of diffusion 
parameters. For illustration, a global boost factor of 2 has also been included in the signal. 
The antiproton flux is compared to several measurements, whereas the expected statistical 
error after 3 years of data sampling by AMS-02 is indicated. This plot has been borrowed 
from [13]. 

Mate predominantly into W + W~ (79.9%) and ZZ (20.1%) gauge boson pairs. 
For illustration, a global boost factor of 2 has also been included in the DM 
signal. The three last PAMELA data (red) points lie just between the astro- 
physical background (yellow band) and the predicted signal (red solid). The 
error bars are sufficiently large to still allow the presence of a 1.7 TeV Wino-like 
species in the Galactic halo, although the pure astrophysical explanation is per- 
fectly compatible with the measurements. The AMS-02 experiment is now on 
board the ISS, and its large acceptance and impressive statistics should allow 
to disentangle both possibilities, as indicated by the theoretical dark-green data 
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points. 



4. The PAMELA positron excess 

The confirmation by the PAMELA collaboration [U] of a positron excess 
above 10 GeV has been triggering a lot of excitement in the field of particle as- 
trophysics since its announcement three years ago. This excess has been readily 
considered as the first long waited hint of the presence of WIMPs in our Galaxy. 
However, the DM candidates which can potentially lead to this positron anomaly 
must have quite special properties. They are also severely constrained by radio 
and gamma-ray observations, unless they are tightly packed inside improba- 
ble or bizarre clumps. These species could also be unstable with abnormally 
long lifetimes. In order to reach any conclusion, we need to investigate how 
positrons (and electrons) propagate throughout the DH and we must calculate 
the astrophysical background to the PAMELA signal. Although this positron 
excess could be generated by annihilating or decaying DM particles, William of 
Ockham would warn us that a more natural explanation is to be found in pulsars 
for instance, and that entia non sunt multiplicanda praeter necessitatem. 

Jf..l. Positron propagation and associated Green function 

The propagation of CR positrons in the magnetic halo of the Galaxy dif- 
fers from that of nuclei in several respects. Although space diffusion is still 
an essential ingredient common to all CR species, positrons undergo mostly in- 
verse Compton and synchrotron energy losses, as discussed in [45] for instance, 
whereas nuclei and antinuclei are mostly sensitive to the Galactic wind and to 
nuclear interactions as they cross the Milky Way disk. As a result, a positron 
line injected by a source leads to an extended positron spectrum once propa- 
gated. This is at variance with most of nuclear species for which energy losses 
can be neglected in a first approximation. Consequently, the diffusion equation 
that relates the positron space and energy density ip to the source term q(x, E) 
takes the form 

- K (^M + A { b ^s {E) ^ = q{X:E) . (22) 

The first term is simply the diffusion coefficient written as K(E) ss Ko(E j Eq) s , 
where Eo = 1 GeV will be used hereafter to keep track of the correct units. The 
synchrotron and inverse Compton processes lead to energy losses with the rate 

6 ioss (£;) ^ = _^o e 2 (23) 
dt v ' 

where e = E/Eq is the reduced positron energy. The energy loss timescale te 
is set equal to 0.6 x 10 16 s. The master equation (22) can be solved with Baltz 



& Edsjo [46] clever trick which consists in translating the energy E into the 
pseudo-time 



t(E) = T E x ] v(E) = — } . (24) 
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In this formalism, the energy losses which positrons experience boil down to a 



mere evolution in the pseudo-time t and the diffusion equation ( 22 1 simplifies 
into the well-known heat equation 

^ - K Ai> = qf(x,t) . (25) 
at 

The positron density is now tp = e 2 tp whereas the positron production rate 
has become q e + — e 2 ~ s q e + . Both ip and q have the same dimensions as be- 
fore, because e is a dimensionless variable. The Green function G of the heat 



equation ( 25 1 leads to the positron propagator through 



G e+ (x, E <- x s , E S ) = t?^ G (x, t <- x s , i s ) , (26) 
where the connection between the energy E and pseudo-time t is given by re- 



lation (24 1. We are led to express the density of positrons resulting from their 



transport within the Milky Way as the convolution 

rE s = + oo p 

ip e+ {*,E) = dE s d 3 x s G e+ (x,E ^x s ,E s ) q e+ (x s ,E s ) , 

JE S =E iDH 

(27) 

recovering thus the generic expression ([8| . The positron propagator G e +(x, E <— Xs 
measures the probability for a positron injected at xg with energy Eg to reach 
the location x with the degraded energy E < Eg. In the 3D limiting case of an 
infinite magnetic halo, the heat Green function connecting the source xs to the 
Earth is the Gaussian distribution 

G(x ,^x s J s ) = 9(f) {^pxP^^} , (28) 

where f = t — t$ is the typical duration, including the diffusion process, over 
which the positron energy decreases from E$ to E. The distance between the 
Earth and the source is denoted by 

r s = {(x Q -x s ) 2 + (y e -ys) 2 + (ze-z s ) 2 } 1/2 . (29) 



The radial behavior of the heat propagator ( 28 1 suggests to define the character- 
istic diffusion length Ad = V^K f which sets the scale of the positron sphere, 
i.e., the region where most of the positrons detected at the Earth are pro- 
duced. It depends on the injected E$ and detected E positron energies through 
the pseudo-time difference f. Above GeV energies, Ad is typically smaller than 
5 kpc. Most of the positrons detected by high-energy experiments like PAMELA 
have a local origin. 

The magnetic halo is nevertheless finite, with radial and vertical boundaries. 
Taking into account the former requires to use an expansion over Bessel func- 
tions as explained in [47]. However, in most situations, it is safe to ignore the 
radial boundaries and to picture the DH as an infinite slab with half-thickness 
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L. Sources located beyond R — 20 kpc should be disregarded though, since 



the convolution (27) is only performed over the DH. The infinite slab hypothe- 
sis allows the radial and vertical directions to be disentangled and the reduced 
propagator G may be expressed as 

G (x, i «- x s , i s ) = . ^ - exp (- } V (z, t <- z s , is) , (30) 

where f = i — is as before. The radial distance between the source xg and 
the point x of observation is now defined as r = {{xq — xs) 2 + (y© — 2/s) 2 } 1 ^ 2 - 
The vertical boundary conditions need to be implemented. Wherever the source 
inside the slab, the positron density vanishes at z = ±L. 

(i) A first approach relies on the method of the so-called electrical images and 
has been discussed in 0S]. Any point-like source inside the slab is associated 
to the infinite series of its multiple images through the boundaries, at z = ±L, 
which act as mirrors. The n-th image is located at z n = 2nL + (—l) n zs, and 
has a positive or negative contribution depending on whether n is even or odd. 
When the diffusion time f is small, the vertical function V is well approximated 
by its infinite ID limit 

V (z, i «- z s ,ts) - Vid (z, t <- z s , t s ) = 6 {T l _ exp J - ^ Z ~ Zs ) \ ■ 

(31) 

A relevant parameter is the ratio £ = L 2 / (4Kqt) = L 2 /Aq. In the regime where 
it is much larger than 1, the propagation is insensitive to the vertical boundaries. 
On the contrary, when £ is much smaller than 1, a large number of images needs 
to be taken into account in the sum 



E 

n — — oo 



V(z,ii-z s ,is) = V\ (-1)" Vid (z,t ^ z n ,t s ) , (32) 



and convergence may be a problem. 

(ii) It is fortunate that a quite different approach is possible in that case. The 



ID version of equation ( 25 1 actually looks like the Schrodinger equation (though 
in imaginary time) that accounts for the behavior of a particle inside an infinitely 
deep ID potential well extending from z = —L to z = +L. The eigenfunctions 
of the associated Hamiltonian are given by 

(fln(z) = sin{fe n (L — \z\)} (even) and (fi' n (z) = sin {k' n (L — z)} (odd) . 

(33) 

They depend on the vertical coordinate z through the wave- vectors k n = (n — 1 /2) 
(even) and k' n = nn/L (odd) . The vertical propagator may be expanded as the 
series 

V(z,i^Z S ,is) = 1 { e ~ A " f <Pn(z S ) Vn{z) + e~ X n f ip' n (z S ) ip'Jz)} , 

n= 1 

(34) 
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where the time constants A„ and A^ are respectively equal to K Q k\ and K k' n . 
In the regime where £ is much smaller than 1, i.e., for very large values of the 
diffusion time t, just a few eigenfunctions need to be considered in order for the 



sum (34 1 to converge. Notice finally that the energies E and Eg always come 



into play in the reduced propagator G through the diffusion length Ar> 

lf.,2. The background of secondary positrons 

Like for antiprotons, an irreducible background of secondary positrons is 
produced by primary CR nuclei colliding on the ISM. The dominant mechanism 
is the collision of protons with hydrogen atoms at rest, producing charged pions 
tt^ which decay into muons ^ . These are also unstable and eventually lead to 
electrons and positrons through the chain 

P + H > X + 7T ± 

M ± — > + v e + e ± . (35) 

In proton-proton collisions, pions can be produced in two different ways, de- 
pending on the energy E p of the incoming proton. Below ~ 3 GeV, one of the 
protons is predominantly excited to a A + resonance which subsequently decays 
into a nucleon and a pion. As isospin is conserved, the branching ratios into 
p + 7T° and n + n + are respectively equal to 2/3 and 1/3. Above ~ 7 GeV, the 



pion production (35) is well described in the framework of the scaling model. 
Various parameterizations are given in the literature [351 E2] f° r the Lorentz 
invariant pion production cross section E v (d 3 a/d 3 p 7r ). Positrons may also be 
produced through kaons generated in proton-proton collisions. In a first chain 



of reactions, the kaon plays the same role as the pion in the set of decays (35). 
The branching ratio of the K ± — > + fjr 1 decay channel is 63.5%. The kaon 
may also decay into a pair of pions, through the reaction K — > 7r° + w , 
with a branching ratio of 21%. The charged pion tt + subsequently follows the 



chain (35) and yields a positron. These two main kaon decay modes contribute 
together a few percent to the total positron production differential cross sec- 
tion. Parameterizations of kaon production in proton-proton interactions can 
also be found in [T^l EE] for the scaling regime. Useful parametric expressions 
for the yield and spectra of stable secondary species produced in proton-proton 
collisions have been derived from experimental data and summarized in |49j . 
Primary CR protons induce a production of positrons whose rate, per unit of 
volume and energy, is given by the integral over proton energy 

f d 
ql% c (x, E e ) = 4 TT n H (x) / $ p (x, E p ) x dE p x ^(E p -> E e ) . (36) 

As in the case of antiprotons, this relation can be generalized in order to incor- 
porate CR alpha particles as well as interstellar helium. The positron produc- 
tion rate g s + c (x, E e ) can be approximated by its solar value since most of the 
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positrons detected at the Earth originate from the solar neighborhood. Folding 



that rate in relation (27) yields the flux 



a r+oc 

$ s e c + c (0, e = E e /E ) - 9^ x 4 x / de s x I (A D ) x <£f(©, e s ) . (37) 



The integral / is the convolution of the reduced positron Green function over 
the Galactic disk alone 

/(A D ) = / d 3 x s G(x 4-x s ;A D ) . (38) 

J disk 

As regards observations, the ATIC [55], HESS [53155] and Fermi [53] collab- 
orations have reported measurements of the total CR electron and positron flux 
at high energy. The Fermi-LAT instrument, in particular, has collected high 
precision data between 20 GeV and 1 TeV, which can be fitted by the simple 
power-law spectrum [5 9) 

$ e± = {175.40 ±6.09 GeV^m-^sr- 1 } £ -(3-°45±o.oo8) , (39 ) 

where e = E/Eq is the reduced lepton energy. Primary CR electrons are believed 
to originate from the interstellar medium, like most of the primary CR nuclei, 
and are accelerated by supernova driven shock waves which inject them with 
a spectral index a ~ 2.2 ± 0.1 inside the Milky Way disk. They subsequently 
propagate within the magnetic halo and lose energy. A detailed investigation [5U] 
leads to a CR electron flux at the Earth of the form 

$ e - oc e~ a - °' 5 - 6 1 2 , (40) 

with a spectral index of order 3.05 ± 0.15, in reasonable agreement with the 
Fermi-LAT measurements should 6 = 0.7, as indicated by the MED model of 
table [I] Some authors [61] find that value too high. It suffices then to increase 
a in order to match the Fermi-LAT observations. The lepton spectrum contains 
also a small admixture of secondary electrons and positrons, produced by CR 
primary nuclei impinging on interstellar gas. These spallation reactions generate 
many charged pions which eventually decay into positrons and, in a lesser extent, 
into electrons. The source spectral index is set by the CR proton and helium 
fluxes. With a gross value of a — 2.7, we expect a spectral index of 3.55 ± 0.05 
for secondary positrons. 

The PAMELA collaboration has measured [H] the positron fraction e + /(e~ + e + ) 
over a large energy range. From the above mentioned arguments, we would 
naively expect that fraction to decrease with energy as e -1 / 2 . A more refined 
analysis |55j confirms that trend, as featured by the banana shape region of 
Fig. [5] The thickness of the yellow band gauges the uncertainties arising from 
CR transport. In spite of these, the positron fraction is expected to decrease 
above 10 GeV. The PAMELA measurements establish on the contrary that the 
positron fraction increases at high energy. This very important observation is in 
contradiction with a pure secondary origin for the positrons. Primary sources, 
which directly inject these particles in the interstellar medium, are necessary. 
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Figure 5: Positron fraction as a function of positron energy. The yellow banana shape region 
encompasses the astrophysical background of secondary positrons yielded by the 1,600 different 
configurations of CR propagation parameters, shown in [9] to be compatible with the B/C 
data. This band is bounded by the MAX (short dashed) and MIN (solid) curves, while the 
central long-dashed curve stands for the MED model (see table 111. The nuclear cross sections 
have been parameterized according to [49| whereas the recent PAMELA measurements of the 
CR proton and helium spectra [50] have been fitted by [51J . In the same figure, the positron 
fraction obtained with the positron flux calculated by 1451 and fitted by 46 a is also indicated by 
the solid black line labeled MS98. The electron and positron flux enters in the denominator 
of the positron fraction. Measurements from AMS-01 1521 and Fermi-LAT |53| have been 
parameterized according to 1541 . A solar modulation with a Fisk potential of 600 MV has 
been applied to the positron flux. This corresponds to the level of solar activity during the 
data taking of AMS-01. This figure has been borrowed from |55| . 



4-. 3. The PAMELA excess as a signature for Galactic WIMPs 

The announcement of a positron excess by the PAMELA collaboration has 
triggered a lot of excitement. This excess was actually considered as the first 
hint of the presence of DM species in the Milky Way halo. Should WIMPs exist, 
their annihilations would be a source of primary positrons. The production rate 
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%+ i given by relation (|2|), yields the positron flux at the Earth 

rm x /E 

$° M (0,e) =Fx^xJ de s f(e 8 ) /dm (Ad) • (41) 

The halo integral /dm |47j is the convolution of the reduced positron propagator 
G with the square of the Galactic DM density 

/dm (A d ) = / d 3 x s G (x Q <- x s ; A D ) (^1) . (42) 
Jdh I P© J 

The positron excess lies above 10 GeV. It can be generated by a particle with 
a mass in the range 100 GeV to a few TeV, in good agreement with theoretical 
expectations for WIMP masses. If the DM species are thermally produced 
during the big bang, their relic abundance matches the WMAP value [63] of 
n x h 2 = 0.1126 ± 0.0036 provided that their annihilation cross section, at the 
time of decoupling, is equal to (<r ann v) ~ 3 x 10~ 26 cm 3 s _1 . The energy 
distribution /(e) of the positrons produced in a single annihilation event is 
set by the WIMP type. High energy positrons (and electrons) cannot diffuse 
on long distances, and those detected at the Earth must have been produced 
locally, hence a DM density of p& — 0.3 GeV cm~ 3 . Baring in mind these 
benchmark values, one can estimate q ^ 1 as well as the positron flux <f> D _^. 
Alas, the signal derived for a large choice of DM candidates, and for various 
DM Galactic distributions, is way too small to account for the observed excess. 
For a WIMP mass m x of 1 TeV, g D + M needs to be enhanced by a factor of ~ 10 3 
in order to match the PAMELA measurements, as illustrated in the left panel 
of Fig. § 

To save the WIMP explanation of the PAMELA positron anomaly, a first 
remedy is the possibility that (<7 ann f ) is much larger than what is currently 
assumed. Three directions at least have been explored so far. To commence, 
WIMP decoupling in the early universe takes place conventionally during a pe- 
riod of radiation domination. If the universe at that time is dominated by 
another component, like a scalar field rolling down its potential [BH [55], the 
expansion rate is significantly increased and the same WIMP relic abundance 
Q x requires now a much larger annihilation cross section. Another possibility is 
that WIMPs do not decouple from the primordial plasma. In this class of less 
natural scenarios, a non-thermal production takes place through the decay of 
heavier species, like gravitinos or moduli fields [BB]. The WIMP relic abundance 
is no longer determined solely by (cr ann i>). It also depends on the production 
rate and lifetime of the unstable heavier states. A final solution relies on the 
so-called Sommcrfcld effect [12 03] . If the DM species are massive enough and 
interact through the exchange of a light particle <j>, they are attracted toward 
each other at sufficiently low velocities. Their wave functions get focalized at 
the interaction point so that the annihilation cross section can be significantly 
enhanced. Scenarios based on that effect have been proposed [571 EH1 EH] , reviv- 
ing the idea of secluded dark matter [70]. WIMPs are still thermally produced 
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in the early universe. Their annihilation cross section at that time is equal 
to the canonical value of 3 x 10~ 26 cm 3 s _1 . As the universe expands, the 
average WIMP velocity /? decreases. The Sommcrfcld enhancement sets in as 
soon as the WIMP kinetic energy is dominated by the interaction energy, i.e., 
when (3 no longer exceeds the effective WIMP coupling constant a'. The cross 
section increases and may be today large enough to account for the PAMELA 
positron excess. The enhancement saturates when the WIMP de Broglie wave- 
length ~ l/{(3m x ) exceeds the range ~ 1/m^ of the interaction. This occurs 
for /3 < m ( p/m x . The annihilation cross section is also resonantly increased 
whenever a' ~ a m x /n 2 , where n is an integer [71] , 

Irrespective of the mechanism enhancing (<7 ann v), WIMPs should not over- 
produce antiprotons [7U |H] . The annihilation channels leading to antiprotons 
must be significantly suppressed. The antiproton to proton ratio measured by 
the PAMELA collaboration [73] does not show actually any excess and is con- 
sistent with a pure secondary origin, as featured in the right panel of Fig. [6j 
Therefore, besides an abnormally large annihilation cross section, DM species 
preferentially annihilate into charged leptons, a feature which is unusual in 
supersymmetry, but more typical of Kaluza-Klein theories. From a phenomeno- 
logical point of view, WIMP annihilations must proceed either directly to lepton 
pairs or, in some models, through the production of the above mentioned light 
mediator bosons <j) which subsequently decay into leptons via 

x + x — > 4> + 4> — ■» i + r i + r • (43) 

The possibility to explain the PAMELA positron excess by leptophilic DM 
particles with enhanced annihilation cross section has been extensively inves- 
tigated. Depending on its magnitude, the injection of high-energy positrons 
and electrons is associated to clear astrophysical signatures [7i] . As they spiral 
along the Galactic magnetic field lines, these particles produce a radio emission 
through synchrotron radiation. They also inverse Compton scatter (ICS) pho- 
tons from the CMB and stellar light. Finally, the charged leptons released by the 
WIMP annihilations can produce final state radiation photons. After an exten- 
sive scan [751 [TBI [77] of the various possibilities, the region of the (m x , (cr a nnw)) 
plane compatible with the PAMELA and Fermi-LAT measurements is found 
to be excluded for most of the leptophilic DM candidates. The most stringent 
limit arises from the synchrotron radio emission produced at the Galactic center, 
a region where the DM distribution may be cuspy, and where magnetic fields 
should be strong. That limit weakens considerably, though, as soon as a cored 
isothermal DM distribution is preferred to the usual Einasto or NFW profile. In 
that case, leptophilic WIMPs are not excluded by radio observations. A more 
reliable tool is provided by the ICS of the positrons and electrons produced by 
WIMP annihilation at high Galactic latitude, a region which is less subject to 
uncertainties than the center of the Milky Way. A careful investigation [75], 
performed only in the case where WIMPs annihilate directly into lepton pairs, 
excludes values of m x larger than 1 TeV. Finally, another model independent 
analysis implies the injection of WIMP annihilation products in the intergalac- 
tic medium at redshift ~ 1000. The recombination process is affected and an 
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imprint is left on the CMB, which may be searched for in the high precision 
WMAP data [73]. Most of the PAMELA compatible leptophilic DM candidates 
still survive that test [80], but will be probed by the Planck mission. 

Another way to increase the production rate of DM positrons relies on the 
existence of substructures inside the smooth DM Galactic distribution. Be- 
cause (p 2 (x)) is always larger than (p(x)) 2 , inhomogeneities tend to effectively 
enhance q 1 ^ 1 , a quantity proportional to the square of the DM density p(x). 
How significant is that enhancement depends on the position of the clumps 
with respect to the smooth component of the DM halo. The inner structure of 
clumps, as well as their mass distribution, are also key factors. Finally, the lower 
mass cut-off and the concentration-to-mass relation of the substructures need 
to be considered. A statistical analysis is mandatory since we live inside one 
particular realization of the clump distribution, to be taken out of an infinite 
set of similar realizations. Using the tools specifically forged for computing the 
odds of that Galactic lottery [5T] , a comprehensive analysis [53] indicates that 
the boost factor to be applied to in the case of a ACDM universe cannot 
exceed at most a factor of 20 at high energy. The correct value is presumably 
much smaller, i.e., of order unity. 

Although a clumpy DM halo does not lead, on average, to an increase of 
the positron flux at the Earth, the statistical variance of the boost factor is 
large at high energy, hence the existence of configurations where q ^ 1 could be 
significantly enhanced. The possibility that a DM clump lies, for instance, in 
the vicinity of the solar system has been suggested [531 El] as an explanation 
for the PAMELA positron anomaly. Such a local substructure would actually 
outshine the rest of the Galactic DM distribution. Alas, the probability that 
such a situation occurs in the DM halo of the Milky Way is vanishingly small 
as demonstrated in an analysis [55] based on the results of the cosmological 
N-body simulation Via Lactea If 86J. In the favorable case of a 100 GeV DM 
species annihilating into pairs, the PAMELA positron excess is best fitted 
by a subhalo with annihilation volume 

i= f J^Ml 2 d 3 x _ 1U kpc 3 ( (44) 

J clump I P<3 J 

located at 1.22 kpc from the Earth. That configuration has a probability of only 
0.37% to occur. Other arrangements are even less probable. 

Finally, the possibility of unstable DM candidates with very large lifetimes 
has also been explored (see for instance [57] [55, 89J). Irrespective of the mech- 
anism which produces them in the early universe, unstable species could in 
principle generate the PAMELA signal. The production rate |2]) can now be 
written as 

tff w = {^} M • («> 

where the quantity r\ (a ann v) (p/m x ) has been replaced by the decay rate t7 . 
For a 1 TeV WIMP, a lifetime of order 2 x 10 26 sec is required to reproduce the 
positron excess. The synchrotron radio emission from the Galactic center is no 
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longer a problem, even in the case of a NFW or Einasto profile, since the positron 
production rate from decaying DM is proportional to the DM density p, and 
not to its square. A detailed analysis [78] excludes however unstable DM as it 
would unacceptably contribute to the isotropic gamma-ray emission inferred 
from the Fermi data. This limit does not apply to WIMPs predominantly 
decaying into muon pairs. The lifetime needs to be fine tuned though. Higher 
dimension operators must be invoked to explain the large value of Td ec required 
by PAMELA. If a four-fermion point-like coupling is responsible for the WIMP 
instability, we anticipate [90] a decay rate cx m^./M 4 , where M denotes the 
typical scale of the underlying high-energy theory. If factors of 2 and tt are 
neglected, a GUT mass M ~ 2.3 x 10 16 GeV could lead to the required lifetime. 
Notice that the DM species must also be leptophilic to prevent antiprotons from 
being overproduced. 

The CR lepton anomalies reported three years ago, and the PAMELA positron 
excess in particular, have triggered a boiling activity in the field of particle astro- 
physics. The possibility that these signals are produced by WIMPs annihilating 
or decaying in the Milky Way halo is fading away. The DM explanation is fairly 
contrived and far from the main stream. The positrons which the DM species 
are supposed to produce generate in turn radiations which are not seen. As 
William of Ockham would have pointed out, a more plausible explanation has 
to be found elsewhere. Actually, a simple one lies in pulsars. These objects are 
known to exist, and they do inject (primary) positrons in the Galactic disk. The 
lepton anomalies can be nicely explained by local pulsars and supernova rem- 
nants, as pointed out by [91] as well as by a recent and very comprehensive [60] 
study. 

5. Antideuterons or the new challenge 

Antideuterons are the nuclei of antideuterium, the element that mirrors deu- 
terium. Although they have been seen in nuclear or lepton interactions at ter- 
restrial accelerators, antideuterons have not yet been detected in the cosmic 
radiation. A small admixture is nevertheless expected to be astrophysically 
produced by the spallations of CR protons and helium nuclei on the ISM. Any 
process that generates at the same time an antiproton and an antineutron could 
lead to the production of antideuterons. The requisite is the fusion of the an- 
tiproton and antineutron pair into an antideuteron. This happens whenever the 
momenta of the antinucleons are so close to each other that the kinetic energy 
of the pair, as seen in its rest frame, is of the same order of magnitude as the 
antideuteron binding energy. 

To quantify this effect, let us consider the collisions of CR high-energy pro- 
tons on hydrogen atoms at rest. The number d/Vx of antiprotons, antineutrons 
or antideuterons generated in a single reaction, with momenta kx, is related to 
the differential production cross section through 

dN x = — dV x (Vi,k x ) , (46) 

Otot 
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where otot denotes the total proton-proton cross section. The total available 
energy is ^/s. The corresponding differential probability for the production of 
the species X is defined as 

dAf x = Jx(^i,kx)d 3 k x ■ (47) 

For antiprotons, for instance, it can be expressed in terms of the Lorentz invari- 
ant cross section through 



d 3 Op 

otot Ep Jp(vs,kp) = Ep 



(48) 
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That cross section is experimentally well-known and a parameterization has been 
given, for instance, in |14j . Isospin symmetry can be assumed, furthermore, to 
hold. If so, the production of antineutrons is the same as for antiprotons. The 
calculation of the probability for the formation of an antideuteron can now 
proceed in two steps. We first need to estimate the probability for the creation 
of an antiproton-antineutron pair. Then, the antinucleons merge together to 
yield an antinucleus of deuterium. As explained in |92] , the production of two 
antinucleons is assumed to be proportional to the square of the production of 
one of them. The hypothesis that factorization of the probabilities holds is 
fairly well established at high energies. For spallation reactions, however, the 
bulk of the antiproton production takes place for an energy ^/s ~ 10 GeV which 
turns out to be of the same order of magnitude as the antideuteron mass. Pure 
factorization should break in that case as a result of energy conservation. It 
needs to be slightly adjusted. Following [HUES], we ma y assume that the center 
of mass energy available for the production of the second antinucleon is reduced 
by twice the energy carried away by the first antinucleon 

^p\n(Vs,kp,k n ) = - .Fp(yi,kp) J- n (7i-2Sp,k n ) + (p«2) . (49) 

Once the antiproton and the antineutron are formed, they combine together to 
give an antideuteron with probability 

j- S (^,k S )d 3 k Q = y'd 3 kpd 3 k n c(kp,k n ) j- P!ll (^,k p ,k n ) . (50) 

The summation is performed on the antinucleon configurations for which kp + 
k n = kg. The coalescence function C(kp,kn) describes the probability for a 
p-fi pair to yield by fusion an antideuteron. That function depends actually 
on the difference kp — k n = 2 A between the antinucleon momenta so that 



relation (50 1 becomes 



Jb(^,k ) = |d 3 AC(A) J-p, n (v / i,k p = 1 ^+A,k n = 1 J-A) . (51) 

In the center of mass frame of the proton-proton collision, an energy of ~ 3.7 
GeV is required to form an antideuteron, to be compared to a binding energy 
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of B ~ 2.2 MeV. The coalescence function is therefore strongly peaked around 
A = and expression (51 1 simplifies into 



J- s (^,k 5 )^||d 3 AC(A)|xJ-p ! n(^,kp = ^,k fi = ^) , (52) 

where the probability for the formation of the p-n pair has been factored out. 
The term in brackets may be estimated in the rest frame of the antideuteron 
through the Lorentz invariant term 

%-*AC{A)*[-^-\ . CM) 
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In that frame, as already mentioned, the antinucleons merge together if the 
momentum of the corresponding two-body reduced system is less than some 
critical value P C oai- That coalescence momentum is the only free parameter of 
this simplistic factorization and coalescence scheme. Theoretical values range 
from y/nipB ~ 46 MeV, naively derived from the antideuteron binding en- 
ergy, up to 180 MeV as would follow from a Hulthen parameterization of the 
deuterium wave function (33]. We therefore expect F coa i to lie somewhere in 
the range between 50 and 200 MeV. Direct comparison with accelerator data 
yields a phenomenological value of 58 MeV as shown in [p_2 [553] . A recent up- 
grade [HI [19] points toward 79 MeV. The Lorentz invariant cross section for 
the production of antideuterons can now be derived under the form 

= (t^ps ^ 1 

7-i 



D d3 kB [mpm^J \3 coa V " 2a t . 



x |Sp^(^,kp) x£; fi ^(v^-2^,k fi ) + (pon)j . 

(54) 

A summation over the production angle 9 in the Galactic frame leads to the 
differential cross section 



Hit- f 6m rPr 

^{p(£ p ) + H — s- D(£ s )} = 2 tt fc 5 J E 



d(-cos6) . (55) 
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In that frame, 9 denotes the angle between the momenta of the incident proton 
and produced antideuteron. It is integrated up to a maximal value 9m set by the 
requirement that, in the center of mass frame of the reaction, the antideuteron 
energy E^ cannot exceed the bound 

s — 16 + mi 
Eh M = • (56) 



The integral (55) is performed at fixed antideuteron energy E^ — mj^ + fcg. 
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Figure 7: The various contributions to the antideuteron flux at the Earth, modulated at 
solar minimum, are plotted as a function of kinetic energy per nucleon Tk. The astrophysical 
background (black dashed line) is generated by the spallations of CR protons and helium 
nuclei on the ISM. This secondary component is featured for the MED model of table ^ 
Below a few GeV/n, the interactions of CR antiprotons with the ISM need also to be taken 
into account. The DM signal (red solid curves) is derived for a m x = 50 GeV WIMP and for 
the three propagation models of tablejl] The upper dashed horizontal line features the current 
BESS upper limit on the search for CR antideuterons. The three horizontal solid (blue) lines 
are the estimated sensitivities of forthcoming measurements which will be performed by (from 
top to bottom) AMS-02 [52] [95] and GAPS on a long (LDB) and ultra-long (ULDB) duration 
balloon flights [96] [S3 [98] [99] This figure has been borrowed from [19]. 



The above factorization and coalescence scheme can also be applied to WIMP 
annihilation. In that case, the differential multiplicity for antiproton production 
is expressed as 

dNf, „ m dN£ 

4 - £ fl S' 4 • < 57 > 

P F,h P 

The annihilation proceeds, through the various final states F, toward the quark 
or the gluon h with the branching ratio B^ h . Quarks or gluons may be directly 
produced when DM species self annihilate. They may alternatively result from 
the intermediate production of a Higgs or gauge boson as well as of a top quark. 
Each quark or gluon h generates in turn a jet whose subsequent fragmentation 
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and hadronization yields the antiproton energy spectrum dN^/dEp. Because 
DM particles are at rest with respect to each other, the probability to form, say, 
an antiproton with momentum kp is essentially isotropic 

dNf, 

^(X + X^P + ---) = 47rk p EpTp(Vs = 2m x ,Ep) . (58) 

The factorization and coalescence scheme leads to the antideuteron differential 
multiplicity 

<™D _ fiP,L,\ ( ">D \ V- D (F) fdW « 



oLE c V 3fc D / \mpm n j — - ^ u^ p ^ 

It may be expressed as a sum, extending over the various quarks and gluons h, as 
well as over the different annihilation channels F, of the square of the antiproton 
differential multiplicity. That sum is weighted by the relevant branching ratios. 
The antineutron and antiproton differential distributions have been assumed to 
be identical. 

The black dashed line of Fig. [JJ features the astrophysical contribution to 
the antideuteron flux at the Earth. It has been derived with the MED model 
of table [T] Varying the CR propagation parameters does not change much this 
contribution, with a spread in possible values reaching at most ± 40-50 % at 
energies below 1 GeV/n and decreasing down to ~ 15 % above 10 GeV/n, as 
shown in [T(5] . Secondary antideuterons are not much affected by CR propaga- 
tion uncertainties. We already found a similar trend for secondary antiprotons. 
Both species are produced by the interactions of primary CR nuclei on the ISM, 
a process reminiscent of the fragmentation of CR carbon nuclei into boron as 
they collide on interstellar gas. Constraining CR propagation from the B/C ra- 
tio proves to be a powerful mean of predicting the flux of background secondary 
antiprotons and antideuterons. Such is not the case for the DM signal. The 
example of a 50 GeV neutralino is presented in Fig. [JJ for the three CR prop- 
agation configurations of table [T] The spread of the red solid curves reaches 
now two orders of magnitude. The larger the half-thickness L of the DH, the 
stronger the signal. The magnetic halo works like a fisher's net inside which 
the products from DM annihilation are entrapped. As noticed in [U3], the DM 
signal dominates at low energies, below a few GeV/n, over the background. 
This may potentially help to discriminate it from its astrophysical background, 
since very few secondary antideuterons are expected in that energy range. Ac- 
tually, these are produced by energetic CR nuclei impinging upon atoms at rest 
of the ISM. The center of mass frame of the reaction moves with respect to 
the Galaxy. Even though a substantial amount of secondary antideuterons are 
produced with small velocities in this frame, few remain at low energy after 
the proper boost takes them back to the Galactic frame. On the contrary, DM 
species annihilate at rest with respect to the Milky Way, hence a dominant 
contribution of primary antideuterons below a few GeV/n. 

A few processes, which were not considered in the exploratory work of [93], 
have been since then suspected to flatten the spectrum of antideuterons by 
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replenishing its low energy tail, erasing the above mentioned discriminating 
difference between the primary and secondary components. Inelastic but non- 
annihilating (INA) interactions of antideuterons with the ISM may actually 
lead to such a flattening. In principle, these collisions are expected to essentially 
break the antideuteron whose binding energy is very small. That is why the INA 
cross section was disregarded in [53] . But laboratory measurements of the INA 
cross section of antiprotons colliding on deuterons yield [TH] a non-vanishing, 
albeit small, value of order 4 mb at most. Another source of flattening comes 
from CR antiprotons impinging on the ISM and producing antideuterons. 
This process cannot be ignored below ~ 1 GeV/n, as featured in Fig. 1 of [T5]. 
Those two mechanisms (INA interactions and production from CR antiprotons) 
generate tertiary antideuterons and a spectrum at low energy somewhat flatter 
than previously anticipated [93] . The detailed investigation carried out in [19] 
confirms nevertheless the possibility of discriminating the DM signal from the 
astrophysical background, as is clear in Fig. [7J Low energy antideuterons provide 
then a unique tool to probe the presence of DM particles within the Milky Way 
halo. The three horizontal solid blue lines of Fig. [7] indicate the potential reach 
of the forthcoming measurements which will be performed by AMS-02 [52 [55] 
and by GAPS. The later is a detector with a very high rejection rate against 
antiprotons. It will be flown on long (LDB) and ultra-long (ULDB) duration 
balloon flights around Antartica [55 ] [971 19"51 19"9 ] . 

The hypothesis that factorization holds is certainly conservative. We have 
naively assumed that both constituents of the antideuteron are independently 
and isotropically distributed in the center of mass frame of the reaction. This 
is certainly true for the first antinucleon and its associated jet. However, once 
the axis of the pair of jets is determined, the second antinucleon tends also to 
be aligned along that direction. Assuming that spherical symmetry holds in 
that case leads to underestimate the probability of fusion. If both antinucleons 
are back to back, they do not merge actually. But if they belong to the same 
jet, their angular correlation is stronger than what has been assumed above, 
hence an enhanced probability of fusion. The naive factorization and coales- 
cence scheme presented in this section has been recently improved by taking 
into account the collimation of the antinucleons produced inside the same jet. 
The larger the jet energy, the better the collimation and the easier the fusion. 
The momenta of antinucleons appear to be more parallel in the Galactic frame 
than in the jet frame, as a result of the Lorentz transformation that connects 
both systems of coordinates. Modeling hadronization inside jets with PYTHIA 
allows to take more accurately into account the angular correlations between 
the antideuteron constituents than in the naive isotropic scheme. The proba- 
bility of fusion increases [100] typically like the square of the jet energy. The 
secondary component is mostly generated at low y/s and is not affected by this 
effect. On the contrary, the collimation of the antinucleons produced by WIMP 
annihilation leads to a stronger DM signal then previously estimated. The pri- 
mary component no longer decreases like l/m^ and may even overcome the 
background above tens of GeV, should the DM species be very heavy. Notice 
however that a mass in excess of ~ 10 TeV does not seem very realistic. Be- 
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sides, such a DM species would still evade |101j the potential reach of the above 
mentioned forthcoming experiments. 



6. Hot spots in the gamma-ray sky 

6.1. Production mechanisms and DM line of sight integral 

WIMP annihilations also generate high-energy photons whose energy distri- 
bution is described by the function /(-Ely) = dN 1 /dE 1 . The corresponding flux 
at the Earth, from the direction toward which the unit vector u is pointing, is 
given by the product 



x 



los 



This formula is often seen as the emblem of particle astrophysics insofar as it 
exhibits two distinct pieces. The first part is related to particle physics and 
encodes information on the WIMP properties such as its mass and annihilation 
cross section. The second term is clearly astrophysical in nature and deals with 
the distribution of DM along the line of sight (los) toward which u is pointing. 
As regards the particle physics aspect, three different contributions to f(E 7 ) 
need to be considered. 

(i) The dominant source of high-energy photons is related to the production of 
quarks and gauge bosons which subsequently fragment and decay into secondary 
gamma-rays through essentially the two-photon decays of neutral pions 

X + X~> qq,W+W-,...^~f + X . (61) 

For each annihilation channel, this leads to a continuum whose spectrum has 
been parameterized in 102 with the generic form 

dN cout 

— 1 — = exp (a + bx + cx 2 + dx 3 ) , (62) 

dx 

where x = E 1 /m x . This distribution exhibits a characteristic x^ 1 ' 5 power law 
behavior for small values of x and a smooth cut-off when the photon energy is 
close to the WIMP mass. 

(ii) A particularly clear signal of the presence of DM species inside the Milky 
Way halo is the production of monochromatic gamma-rays |103j through the 
reaction 

X + X ^ 7 + 7 & 1 + Z° . (63) 

This process gives rise to characteristic line signals which cannot be mistaken 
for some conventional astrophysical source and which would unequivocally sig- 
nal the presence of an exotic component inside the Galaxy, should a peak be 
detected in the high-energy spectrum. The energy of the photons produced in 



reaction (63) is respectively equal to E 7 = m x and E 7 — m x — (m z /Am x ). 
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Because WIMPs are electrically neutral, the production of monochromatic pho- 
tons is necessarily mediated by loop diagrams. It is generally suppressed and 
the integrated photon yield amounts to ~ 10 -3 of the total. This leads to a 
clear but faint signal which is beyond the reach of current detectors unless the 
process is efficiently enhanced, as in the case of heavy Wino-like neutralinos 
discussed in [251 1351 ITMj . 

(iii) Finally, as already pointed out in |105j . a single photon may be produced 
through internal bremsstrahlung as the WIMP pair annihilates. The gamma-ray 
is radiated by the charged particles that are either exchanged (virtual internal 
bremsstrahlung or VIB) or produced in the final state (final state radiation or 
FSR). This process becomes particularly important for a sizable branching ratio 
into electron-positron pairs, as in the case of MeV DM |106| or Kaluza-Klein 
inspired models [107] . The FSR spectrum associated to the production of the 
charged lepton pair l + l~ is, in leading logarithmic order, well approximated by 
[TU71fT0"H] 



dN* SR = d(a l+l - J v)/dx ^ (x ^ , ^ ^ , ^ 

dx &i+i- v TT ■ 1 




(l-x)\ . (64) 



Final state radiation produces more photons than fragmentation does near the 
upper edge E y = m x of the spectrum and is thus responsible for a characteristic 
sharp cut-off there. For Wino-like heavy neutralinos, photons are radiated by 
final W + W~ pairs as discussed in [109] . 

(iv) Photons can also be produced, albeit indirectly, by the inverse Compton 
(IC) interactions of CR electrons and positrons on the Galactic radiation field. 



The resulting gamma-ray flux is no longer given by relation (60) but depends 
on the densities of CR electrons and positrons. These species will be generically 
designated here as electrons. The IC process is very important for DM species 
which preferentially annihilate into charged leptons, and most of the models 
designed to explain the PAMELA positron anomaly are excluded by the over- 
production of IC photons to which they lead. The scattering of an incoming 
electron with energy E e (and Lorentz factor -f e = E e /m e ) on a radiation field 
of incoming photons with energies E- m generates the gamma-ray spectrum [110] 

3<7t f 1 dq diii 

4 7e 7l/4 7 | 1 dE u 



P i (K v E e ,x)= -£/ ^ -^{E in (q),x} 



*** + < + 1 -* + w^} ■ (65) 

where ctt = 0.665 barn is the Thomson cross section. The parameter q is related 
to the energy E- lri of the radiation field photon through 

q = — E ~< (66) 

r e (E e — E 7 ) 

The dimensionless parameter T e = 4i?i n 7 e /m e determines the regime of the IC 
scattering, i.e., either ultra-relativistic (r e 3> 1) or non-relativistic (r e -C 1) 
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in the Thomson limit. Incoming photons belong to the CMB, the IR radiation 
field and the light from stars. These three contributions are represented by the 



subscript i in relation ( 65 1 . The flux at the Earth of IC photons is given by 
the convolution of the various contributions Vi with the CR electron density 
ip e = dn e /dE e along the line of sight 

1 f f m>c . 
4>\ c (E 7 ,u) = — I ds dE e Xi/j e (E e ,x)x22'Pi( E j> E e,X-) ■ (67) 

J los J m e i 

The CR electron density ip e is generated by the collisions of CR protons and 
helium nuclei on the ISM (background) and also by WIMP annihilations (DM 
signal) . It needs to be calculated everywhere inside the Galactic magnetic halo. 



In the case of the two-photon line, relation (60) leads to the estimate 



$^ M (£L = m Y ,u) = 1.88xl0~ i3 photons cm" 2 s" 1 sr" 1 x ' 'f x J(u) 
7 m wo 2 



The annihilation cross section (cr 11 v) and the WIMP mass m x are respectively 
expressed in units of 10 -29 cm 3 s _1 and 100 GeV. A Majorana type DM species 
has been assumed here with r\ = 1/2. The los integral J(u) has been redefined 
as 

J(u) = {rItq}- 1 x f p 2 (x) ds , (69) 

J los 

with a solar neighborhood DM density of po = 0.3 GeV cm -3 . The galacto- 
centric distance r© of the solar system has been set equal to 8.5 kpc. The los 
integral J depends on the Galactic DM distribution. The various halo models 
of table [2] yield similar values all over the sky, except in the direction of the 
Galactic center where predictions can vary by several orders of magnitude, de- 
pending on the assumed profile. As shown in 111 , where it has been averaged 
over a solid angle of 10 -5 sr, J is respectively equal to 30 and 1.45 x 10 4 for 
the cored isothermal and NFW distributions. In the case of a Moore profile 
which has been adiabatically compressed by the collapse of the supermassive 
black hole lying at the center of the Milky Way, J may even reach the extreme 
value of 3 x 10 8 . The astrophysical uncertainties are enormous. The possibility 
of a strong signal has triggered a febrile activity around the Galactic center 
and motivated observations focusing in that direction. Let us concentrate, for 
instance, on the NFW profile of table [2] For small galactocentric distances, the 
DM distribution simplifies into 

p(r) = H P0 { r f} , (70) 

where the normalization constant H ~ 2.03. The los integral J is found equal to 
7rH 2 /a at an angular distance a from the Galactic center. Once averaged over 
a disk with angular radius 9, this leads to (J) = 27rH 2 /a and to a numerical 
value of - 1500 x (l°/9). 
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6.2. Atmospheric Cherenkov telescopes - a new generation of instruments 

High-energy photons can be detected by satellite borne instruments orbiting 
the Earth and free, therefore, from the screening of the atmosphere. Because 
of the reduced payload which can be carried up in orbit, the collecting area of 
gamma-ray telescopes like Fermi-LAT^] is quite reduced. Above a few tens of 
GeV, atmospheric Cherenkov telescopes (ACT) come into play and offer a nice 
alternative which complements space observations. The effective detecting area 
can be actually quite large insofar as any high-energy gamma-ray impinging 
upon the upper atmosphere generates a shower, and is degraded into many 
optical photons which ordinary telescopes spread on the ground can detect. 
The tracks left in the focal planes of these telescopes give the direction of the 
shower on the sky. Several instruments allow the stereoscopic reconstruction 
of the direction of the initial particle. The Cherenkov light illuminates on the 
ground a disk whose diameter is of order 250 m, hence a large effective area of 
~ 0.05 km 2 for each telescope, to be compared with the typical square meter 
collecting surface of a space instrument. The intensity of the image depends on 
the energy of the incident photon. The ACT technique is thus a powerful tool 
allowing the observation of the high-energy gamma-ray sky. 

However, as for any other observation, a gamma-ray signal is only detectable 
as long as it emerges from the background. More specifically, the signal has to 
be larger than the statistical fluctuations, also called the noise, of the latter. 
For atmospheric Cherenkov telescopes, the dominant background arises from 
CR electrons that penetrate the atmosphere inside which they produce electro- 
magnetic showers of the same type as those induced by high-energy photons. 
It is not possible to distinguish photons from electrons since both species lead 
to the same light pattern on the ground. Fermi-LAT measurements of the CR 



electron and positron flux can be fitted [59] by the power-law ( 39 1 . Imping- 
ing CR hadrons also interact with the atmosphere. The showers which they 
generate tend to develop at a lesser altitude and are more widely spread on 
the ground than those of the electromagnetic type. Stereoscopy is a powerful 
tool to discriminate hadrons from electrons and gamma-rays since the pattern 
recognition of the light pool is then possible. Observations performed between 
50 GeV and 2 TeV yield a hadron flux [TT1] 

$had(£) = 1-8 GeV" 1 cm" 2 s" 1 sr" 1 (E/l GcV)~ 2 ' 75 . (71) 

A small fraction of the hadron-induced showers are mistaken for gamma-ray 
events though. This is the dominant source of background at high energy as 
hadrons have a harder spectrum than electrons. The HESS collaboration quotes 
a rejection factor of one misidentified event over a sample of 300 showers gen- 
erated by CR protons. Satellite borne instruments do not suffer from the same 
problem. However, point sources can still be buried inside a Galactic gamma- 
ray diffuse emission produced, above 100 MeV, by the spallation of interstellar 



3. See the web site http://www-glast.stanford.edu/ and references therein. 
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gas by CR nuclei. The flux of this diffuse emission is given by the convolution 
along the los of the hydrogen density iih with the gamma-ray emissivity In per 
hydrogen atom 

$ s 7 oc (£ 7 ,u) = f n H (x)xJ H (x,£ 7 )xrfs . (72) 

J los 

The Galactic diffuse emission is made of secondary photons resulting mostly 
from the interactions of high-energy CR protons and helium nuclei on the ISM. 
This process has already been discussed for antiprotons and positrons. It has 



been respectively described by the source terms q^ c (x,Ep) in relation (12) of 



section 3.1 and (?g+ c (x, E e ) in relation (36) of section 4.2 The photon emissivity 
per hydrogen atom is defined as 

g 7 cc (x,£ 7 ) 

2h(x,£ 7 ) = -j -— , (73 

4 7r nn(x) 

and may be written as the convolution over proton energy of the CR proton 
flux with the differential photo-production cross section of proton-proton inter- 
actions 

f do 
/ H (x,£ 7 ) = J $p(x,£ p ) x dE p x j°r(E p -> £ 7 ) . (74) 

The emissivity Ih is expressed in units of GeV -1 s _1 sr _1 since it is essentially 
a production rate per unit of energy and solid angle. In the solar neighborhood, 
the effective gamma-ray emissivity of the ISM, whose atoms are illuminated by 
the local CR protons and helium nuclei, may be approximated, between 3 GeV 
and 1 TeV, by the power law |113j 

Iif(©,£ 7 ) = 3.55 x 1(T 27 GeV- 1 s" 1 sr" 1 (E^/l GeV)" 2 ' 76 . (75) 

This effective emissivity is scaled back to each hydrogen atom and can be readily 



used in expression ( 72 ) to infer the photon flux. The Galactic diffuse emission 
dominates over an extragalactic component which has been measured |114j with 
the EGRET instrument on board the CGRO satellite 

<E> 7 S (£ 7 ) = 7.32±0.34xl0~ 9 MeV" 1 cm~ 2 s" 1 sr" 1 (S 7 /451 MeV) _210±0 03 . 

(76) 

An atmospheric Cherenkov telescope of the HESS caliber has an effective 
detecting area S of order 0.1 km 2 with four mirrors spread on a 300 m x 300 m 
square. One of the main targets of HESS is the Galactic center where the 
putative WIMPs might have collapsed, producing a hot spot in the gamma-ray 
sky. We will assume in what follows that the DM profile is given by the NFW 
distribution of table [2] The effective time T during which the observation of 
the Galactic center is performed, disregarding the periods of daylight as well as 
the nights during which the Moon shines, will be taken to be a month. We infer 
an approximate acceptance of 

S xT~ 0.01 km 2 yr . (77) 
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Let us also assume that this HESS type telescope monitors a circular field of 
view with angular radius 6 ~ 1° surrounding the center of the Milky Way. For 
illustration purposes, the self-annihilation of WIMPs into photon pairs is the 
only process in which we will be interested here. It produces monochromatic 
gamma-rays with energy = m x . The number of line photons collected during 
such a run is 

N™ = 850 photons x x (!) . (78) 

1 m W0 2 {1°) 

These monochromatic gamma-rays are detected within some energy bin whose 
width is set by the resolution of the telescope. We may safely take an energy 
resolution a(E)/E of order 10%, to be compared to the value of 15% in the 
case of HESS. The energy bin that contains the line has thus a width AE 7 of 
order 0.1 x m x . The better the energy resolution, the narrower the line bin and 
the more visible the peak in the photon spectrum. As discussed above, all the 
bins are filled up predominantly by misidentified hadron and electron events. 
The monochromatic signal from annihilating DM species is detectable only if 
it exceeds the statistical fluctuations of that background. If electron induced 
showers are assumed to be the only source of background, although this may 
not be true at high energy where hadrons come also into play, the number of 
background events collected inside the line bin, during the run, amounts to 

N* K = 4.3 x 10 5 photons x (m 100 )" 2 ' M5 x j!J . (79) 

The DM line signal N® M is deeply swamped into the background N® K and 
seems hopelessly out of reach. However, electronic and hadronic events are 
homogeneously spread on the sky since these cosmic radiations are isotropic at 
the Earth. Changing the direction toward which the telescope is pointing does 
not affect the number of background events. On the contrary, the line signal 
disappears as soon as the field of view no longer encompasses the Galactic center. 
By alternatively pointing the telescope on and off the source makes it possible 
to detect the line signal provided that it exceeds the statistical background 
fluctuations. A good estimate for these is given by the Poisson noise (A^) 1 / 2 . 
Detection of a signal with a significance of n is therefore achieved whenever the 
signal to noise ratio is equal to 

N™/^N™ = n . (80) 

The HESS like telescope of this example would detect at the 3 a level a NFW 
distribution of DM species at the Galactic center should their mass and two- 
photon annihilation cross section fulfill the condition 

(ct 77 d) > 2.3 x 10~ 29 cm 3 s" 1 x (mi 00 )°- 98 . (81) 

The HESS collaboration has observed the Galactic center with unprecedented 
accuracy |115| 1116] above 160 GeV. A strong point-like source called HESS 
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J1745-290 is detected at the positions of the supermassive black hole Sagittar- 
ius A* and supernova remnant Sagittarius A East. Its spectral index is ~ 2.25. 
An important gamma-ray diffuse emission is also seen in the same direction, 
with a similar spectral index. It is clearly associated with the band of molec- 
ular clouds lying in the central region and mapped from their CS line. These 
clouds have been recently penetrated by CR protons and nuclei accelerated by 
a nearby supernova event. The correlation between the intensity of the TeV dif- 
fuse emission and the gas column density of the clouds is striking and suggests 
a uniform density of cosmic rays. No line is observed though. Moreover, the 
gamma-ray spectrum is too hard to be compatible with a WIMP annihilation 
signal. As shown in |107j . a Kaluza-Klein DM species could still yield the same 
flat energy distribution if the contributions from internal bremsstrahlung and r 
lcpton decays are included. The price to pay, however, is an unacceptably large 
mass m x of 10 TeV. A contribution from DM annihilation cannot be ruled out 
either, provided that it contributes less than 10% of the signal. 

6.3. Current observations, constraints and prospects 

The launch of the Fermi satellite in 2009 has considerably improved the 
picture in the field of indirect DM searches. The LAT instrument onboard 
the Fermi satellite has allowed for a tremendous recording of the gamma-ray 
sky with unprecedented statistics and resolution, overtopping the amount and 
quality of the data collected by previous space experiments, like EGRET or 
AGILE. Data taking is expected to last until 2014. The effective detection area 
reaches ^10,000 cm 2 over an energy range spreading from ^30 MeV to ^300 
GeV, with an angular resolution ~ 0.1° in a field of view of ~ 1/5 of the sky. 
Such a low threshold together with the capability of observing the full sky in a 
short time make the Fermi-LAT instrument better than current ground-based 
gamma-ray telescopes for DM searches. The latter are actually complementary 
since they provide a better sensitivity for energies larger than 100-200 GeV, 
but only for rather localized sources (ground-based telescopes are optimized for 
a pointing mode). We still note that much larger Cherenkov telescope arrays 
may reach the required sensitivity in the future to start probing DM signals in 
several classes of targets efficiently |117] . Consequently, the following discussion 
will be strongly biased toward the Fermi-LAT data. 

We will discuss the current results and prospects target by target, reviewing 
some predictions and constraints for DM signals from (i) the Galactic center, (ii) 
the diffuse Galactic and extragalactic emissions, and (iii) extragalactic sources 
(dwarf spheroidal galaxies, neighbor galaxies, and galaxy clusters). We empha- 
size that the observational constraints associated with these targets only concern 
the s-wave part of the annihilation cross section, while the p-wave part is also 
often relevant in setting the cosmological DM abundance. Hence, the canonical 
value of (a ann v) = 3 x 10~ 26 cm 3 /s that will often appear in the following has 
to be taken with caution, keeping the previous remark in mind. 
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6.3.1. The Galactic center 

The Milky Way is obviously the biggest nearby DM reservoir where to seek 
DM annihilation traces. Since DM is expected to concentrate at the centers of 
structures, the Galactic center (GC), which is located at ~ 8 kpc from the Sun, 
should be the most luminous known DM source for the Earthian observer. This 
region was actually observed with the EGRET telescope onboard the CGRO 
satellite, though with a poor angular resolution of ~ 1°, which indeed led to 
the detection a gamma-ray signal in the 0.03-10 GeV energy range |118j . In- 
terpretations in terms of DM annihilation were performed for supersymmetric 
candidates {e.g. [US]), but the main difficulty which arose (but was expected) 
was to distinguish between this putative exotic origin and an astrophysical ori- 
gin of the signal. This points toward the challenging character of looking for 
DM annihilation traces toward the GC: this region is dominated by baryons, 
and high-energy astrophysical processes are known to take place (potential sites 
of CR acceleration, pion production through interactions of CRs with the ISM 
gas, inverse Compton from high-energy electrons, etc.). Anyway, the authors of 
Ref. [119] reached the conclusion that DM annihilation could provide a good 
fit to the EGRET data when taking neutralinos in the 50-100 GeV mass range 
(50 GeV being the lower limit they fixed to their scan over parameters), but 
only for central densities larger than usually expected. Adopting an NFW DM 
profile normalized to p Q = 0.3 GeV/cm 3 at the Sun's position (taken at 8.5 
kpc from the GC), they fell short by 1 order of magnitude, which they argue 
could be compensated by a more spiky profile resulting from baryon cooling 
effects known as adiabatic contraction |120L 11211 1122) . We note that the cur- 
rent state-of-the-art numerical works in cosmological structure formation now 
seem to favor a picture where feedback effects dominate over cooling, making 
contraction scenarios much less motivated (e.g. [123 ). 

The GC was further observed, though at higher energies, with the HESS 
ground-based telescope, and a signal was reported in the 0.25-10 TeV energy 
range |115j . This signal is consistent with a point-like source coincident within 
1' of Sgr A* (the data were taken with an angular resolution of ~ 0.1°). The DM 
interpretation is made difficult because of the hard spectrum, well fitted with a 
power-law of spectral index ~ —2.2 up to 10 TeV, which implies very massive 
DM particles if one assumes their annihilation to be the main photon source. 
The exercise was still performed for supersymmetric as well as Kaluza-Klein 
candidates (e.g. |107j and |109j . respectively), where multi-TeV mass particles 
were found able to saturate the signal, though with a loose fit and to the price 
of an arbitrarily large boost factor of ~ 1000. Such an amplification is generally 
difficult to explain in terms of DM distribution, but non-perturbative effects 
were pointed out to significantly increase the annihilation cross section of such 
heavy particles |43J, though by a factor limited to < 100 to simultaneously 
achieve the correct relic density 124 . Nevertheless, as for the EGRET signal, 
astrophysical scenarios were also shown to be plausible (e.g. [1251 1126] ). making 
the DM scenario less attractive. Despite the absence of a robust smoking gun 
for DM annihilation, the HESS data have still been shown very powerful to 
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Figure 8: Gamma-ray excess found toward the 
GC in 128 129 , and tentative interpretation 
in terms of WIMPs annihilating into pairs of 
tau leptons. The plot is taken from Ref. |129l 
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More recently, several studies were conducted relying on the publicly avail- 
able Fermi data. The authors of Refs. |128[ 1129] have found evidences for an 
extended source centered around the GC, the intensity of which is seemingly 
hard to explain in terms of the point source located at the GC, probably also 
responsible for the TeV emission. These authors have attempted a DM inter- 
pretation of the signal, suggesting that the spatial and spectral fits to the data 
improved significantly by adding a ~ 7—10 GeV WIMP annihilating prefer- 
entially into tauons, and distributed according to a rather spiky profile scaling 
like p oc r -13 (normalized to 0.4 GeV/cm 3 locally). They need an annihilation 
cross section of (a ann v) ~ 7 x 10~ 27 cm 3 /s, rather consistent with cosmologi- 
cal requirements. Moreover, they emphasized that the mass range is consistent 
with the recent hints coming from direct detection experiments — their result 
is shown in Fig. [8] Nevertheless, a criticism to this analysis can be found in 
|130j . where no net excess is found with respect to canonical models of diffuse 
emission. More generally, we note that the GC is a region very complicated 
to model, though as fascinating as challenging, which makes DM searches very 
difficult toward this target unless very specific signatures are observed. Indeed, 
finding excesses with respect to a fortiori non-perfect descriptions of the region 
is not that a surprise. This statement is nicely reinforced by the recent finding, 
in the Fermi data, of large scale bubbles extending away from the GC [131] , to 
which astrophysical explanations have been proposed |131U132lll33] — while the 
believer's eyes could see a DM signal, as apparent as in Fig. 7 of |134j . Hence, 
strong improvements in the understanding of CR acceleration and transport, in 
the modeling of the interstellar material, and in the census of the CR sources 
there, are clearly necessary to perform more reliable and detailed analyses of the 
GC and find DM signals in the near future — complementary with other targets 
and messengers. The encouraging point is that GC observations at gamma- ray 
frequencies still do not exclude a DM contribution. Regarding the experimental 
forecasts, a keypoint will be to improve the angular resolution and sensitivity 
in the 1-100 GeV energy range, which might be partly achieved with the new 
or coming generation of Cherenkov telescopes, like HESS-2 and CTA. 

Finally, for completeness, we may also mention the possibility that DM an- 
nihilation could also contribute to the 511 keV emission extended around the 



41 



GC, which was measured to an unprecedented precision with the Integral satel- 
lite [135] . This would imply very light scalar DM candidates in the MeV mass 
range, which is formally not excluded |136j . Again, many classical alternatives 
exist in astrophysics, which makes it difficult to unambiguously discuss this in- 
tense emission in terms of DM. A specific review on this topic can be found in 
Ref. |137j for further details. 

6.3.2. Diffuse Galactic and extra- galactic emissions 

Fermi-LAT observations of the full gamma-ray sky have allowed to scrutinize 
the properties of the Galactic diffuse emission, as well as those of the isotropic 
extra-galactic component (see e.g. |138l 1139"] ). DM annihilation in the Galac- 
tic halo is also expected to contribute to both these diffusion emissions, which 
makes the all-sky data suitable for searches for departure from conventional as- 
trophysical predictions. This again necessitates reliable models of astrophysical 
foregrounds, and significant efforts are currently made to improve them (see 
e.g. |139l 1113] for the Galactic component, and |140j for the extragalactic one). 

To circumvent spurious effects coming from the potential use of an imperfect 
background model, a first strategy consists in looking for signatures which are 
very specific to DM: gamma-ray lines [1031 EH H42j . or sharp spectral features 
stemming from internal bremsstrahlung diagrams [1061 11091 1143] . The former 
signals are loop suppressed, while the latter come with significant amplitudes, 
inducing a sharp spectral bump close to the WIMP mass. 

Gamma-ray line searches in the Fermi data have been performed by several 
groups, e.g. |144[ I145j . No signal was found, and limits have been derived for 
annihilating as well as decaying DM. The results obtained in Ref. [145], valid 
in the 1-300 GeV energy range, are shown in Fig. [9j The authors assumed an 
NFW halo, with a local density of p© = 0.4 GeV/cm 3 , and a scale radius of 
r s = 20 kpc (variations of the halo shape appear as blue bands in the plots) , and 
derived limits for the annihilation into a pair of photons, and for the decay into 
a photon and a neutrino. It is an easy exercise to recast their results for other 
final states (e.g. jX, where X is any massive particle). The obtained limits 
are not of concern for most of supersymmetric models leading to annihilating 
neutralino DM, while they become stringent for more effective models [146] . 
However, they strongly constrain the degree of R-parity breaking in scenarios of 
decaying gravitino DM, as well as on the stau NLSP lifetime, whenever relevant. 
We refer the reader to the original paper |145j for more details. 

Further limits using three-body annihilation channels including virtual in- 
ternal bremsstrahlung (VIB) photons^jvere derived by the authors of Ref. |147j . 
They assumed an NFW halo profile with p = 0.4 GeV/cm 3 and r s — 20 kpc, 
and parameterized the VIB contribution in terms of the mass splitting between 
the mediator and the DM particle fi = (m m /m x ) 2 (larger than 1 in the R-parity 
conserved supersymmetric context), the larger the values of /i the smaller the 
annihilation rates; though values close to 1 may correspond to co-annihilation 
regions. They selected a large region around the GC to maximize the signal- 
to-noise ratio, and performed a detailed analysis of the diffuse emission. While 
they found a 3.1-cr hint for a signal corresponding to a WIMP mass of ~ 150 
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Figure 9: Results taken from Ref. |145| (the 
results obtained in Ref. I144| are shown as red 
dashed curves). Top left: model-independent 
limits on the gamma-ray line flux as a func- 
tion of energy. Top right: limits in terms of 
DM annihilation cross section into a pair of 
photons. Bottom left: limits in terms of DM 
decay rate (decay into a photon and a neu- 
trino). 



GeV, they adopted a conservative approach in deriving limits on the 3-body 
annihilation cross section. Their results are shown in Fig. |10[ where the blue 
horizontal lines are upper bounds corresponding to the requirement that the 
3-body annihilation cross section (an s-wave), when taken alone, is such that 
the relic abundance Q, x h 2 > 0.1; indeed, 2-body contributions are expected, 
which should reduce the relic density significantly (VIB contributions, though 
s- waves, come with a a cm /n suppression factor). The green curves allow com- 
parisons with the limits obtained from dwarf spheroidal galaxies in Ref. |148] , 
and correspond to requiring (cr a nn")/(5 X 10~ 30 cm 3 /s) < 8 n (m x /GeV) 2 /N* ot , 
where N^ ot is the total number of gamma-rays per annihilation (including the 
2-body contribution). We note that the limits obtained for leptonic channels 
are generally better than those derived from dwarf galaxies. 

Finally, the diffuse gamma-ray flux measured by Fermi |138| I149[ 11501 H51| 
139 has also been used to put absolute upper limits on the continuum gamma- 
ray fluxes expected for most of DM candidates. We may illustrate the main re- 
sults, though not exhaustively, by showing the constraints obtained in Ref. [152] . 
where the authors used the estimate of the isotropic diffuse gamma-ray flux per- 
formed by the Fermi collaboration [150] to constrain DM contributions at the 
Galactic scale (one can also use the extragalactic contribution, though less pow- 
erful and more affected by theoretical uncertainties |152U153lll54j ). In addition 
to the DM induced prompt continuum gamma-ray emissions, the authors in- 
cluded the inverse Compton contributions for all leptonic annihilation channels. 



4. As already mentioned in section |6.l| internal bremsstrahlung photons are made of final 
state radiation (FSR) and virtual internal bremsstrahlung (VIB) photons, the latter being 
associated with the radiation of a t-channel charged mediator (e.g. any sfermion in the super- 
symmetric context). These contributions can be separated in a gauge invariant way |143] . 
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Figure 10: Limits obtained on the 3-body an- 
nihilation cross section for leptonic final states 
plus a VIB photon. These limits were derived 
from the Fermi data and the plots taken from 
Ref. I147| . They are compared to complemen- 
tary limits obtained in Ref. 148 from a study 
of dwarf spheroidal galaxies. 
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Their results are reported in Fig. 11 where the bb, t + t~, and annihi- 
lation channels are considered as representative of most of DM models. The 
plots contain complementary limits from studies of dwarf galaxies and of the 
GC [1551 11561 11571 1158j , and also provide forecast limits for a 5-year Fermi-LAT 
sensitivity. We see that the Galactic diffuse emission leads to constraints at 
the limit of the interesting supersymmetric parameter space, assuming canon- 
ical cosmology and particle phenomenology, delineated as in Ref. |117j . On 
the other hand, as illustrated in the plots, such data make it possible to ex- 
clude more contrived configurations, for example those proposed to explain the 
so-called PAMELA excess [1591 1160] . Many similar results have been reached 
independently, that we cannot, unfortunately, exhaustively cite here. Never- 
theless, we emphasize that a reliable and global picture of the diffuse Galactic 
gamma-ray sky will only come out after important refinements in the astrophys- 
ical background modeling are achieved. These are now necessary for seeking DM 
annihilation traces into more detail. 



6.3.3. Extragalactic sources 

DM is also sought in extragalactic objects, like in nearby dwarf galaxies, 
galaxies, and clusters of galaxies. 

The first class of objects, dwarf galaxies, is likely the most promising in terms 
of DM detection because these targets are strongly DM dominated, with no ex- 
pected high-energy astrophysical phenomena (i.e. background free, though this 
has to be weighted by their positions in the sky where the Galactic foreground 
may be strong). Kinematic data are also available which allow to constrain the 
DM distributions in several systems, and it is shown that when using the cur- 
rent angular resolutions of gamma-ray telescopes at the ground or in space (say 
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Figure 11: Plots taken from Ref. [T52l . The 
dashed lines show the complementary lim- 
its derived from Fermi observations of the 
Draco dwarf galaxy in 1 1551 ; the thin dot- 
ted lines show those obtained by the stack- 
ing of dwarf galaxies in 11561 , except for the 
case, where they correspond to the 
radio constraints found in 159 . The triple- 
dot-dashed lines are limits corresponding to 
the HESS observations of the GC [T571H581 . 
The small green rectangle shows the pa- 
rameter space consistent with the contro- 
versial gamma-ray excess found toward the 
GC in [1281 [129]. The solid line shows 
the Fermi limit without the inverse Comp- 
ton contribution from leptonic annihilation 
channels. The red MSSM/MSUGRA re- 
gions were derived in |117l . The shade pink 
and red regions (/x+ fi~ case) are those which 
are consistent with a DM interpretation of 
the positron fraction measured by PAMELA 
and of the electron+positron spectrum mea- 
sured by Fermi [T591IT60] . 



~ 5'), the theoretical errors expected after the angular average of the squared 
DM density in each object, e.g. stemming from the large uncertainties in the 
inner slope, are actually rather small |161[ 1162] . Several observations of the 
so-called canonical dwarf galaxies (less than a ten of objects) were performed 
in the past by ground-based telescopes, but the limited observation times and 
high-energy threshold induced too weak sensitivities to get interesting limits 
(e.g. [1631 11641 1165j ). As for other targets, the Fermi advent has considerably 
improved the detection potential. So far, yet, no gamma-ray signal originating 
from a dwarf galaxy has ever been observed, leading several groups to compete 
to get the best limits (e.g. [T551 UM[ H151 [TS5] ) . In Fig.[l2] we show the results 
obtained for individual dwarf galaxies in [155] (left panel), corresponding to 11 
months of survey mode operation with Fermi, and those obtained by stacking 
all individual contributions in j!48j . corresponding to about 2 years of Fermi 
data. In the left panel, we see that Draco provides the best limit, though it lies 
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Figure 12: Left: indicative limits on the constrained MSSM parameter space from several 
dwarf galaxies, assuming a generic bb spectrum; taken from 155 . Right: limits after stacking 
individual limits, with a self-consistent treatment of statistical and systematic errors (the 
shaded areas show the systematic errors coming from uncertainties in the DM distribution 
inside the objects); taken from | 148| . 

at about one order of magnitude away from the interesting parameter space. On 
the other hand, the stacking analysis is shown a very powerful alternative, which 
currently yields the best limits on the s-wave part of the WIMP annihilation 
cross section. We notice that despite the rather large systematic errors coming 
from uncertainties in the DM distribution inside the objects (encoded as the 
shaded areas in the plot), WIMP masses below ~ 50 GeV annihilating prefer- 
entially into tau leptons or b quarks at the canonical rate start to be in serious 
tension with the data. Beside exclusion prospects, this is also very encouraging 
in terms of detection prospects. 

As regards expectations for the future, we refer the reader to Ref. |161j . where 
the authors show that future ground-based observatories, like CTA, will allow 
for interesting complementary surveys with respect to Fermi, though further 
observations with Fermi are likely the best strategy for discovery purposes — 
10 years of Fermi data would really drive us in the ballpark, though such a 
duration currently appears as unrealistic in terms of funding (the official end of 
Fermi operations is June 2012, though 2 supplementary years seem granted). 

For completeness, it is important to mention that other extragalactic ob- 
jects might also be interesting sources for DM searches. Neighbor galaxies, 
like M31 (Andromeda), are potentially good targets (1671 11681 I102j . but as 
for the Galaxy, it is rather difficult to predict the astrophysical background in 
these objects, where, at variance with dwarf galaxies, high-energy phenomena 
are expected. Several unsuccessful observations have been performed in the 
past, e.g. [169] 1170] I171j . but the constraints obtained on DM models were in- 
significant. Nevertheless, some giant galaxies, like M87 and NGC 1275 were 
recently observed at TeV gamma-ray energies, but with strong variabilities 
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(e.g. |172l 11731 1174| V The latter objects usually belong to the broader class 
of radio emitters, which allow for rich multiwalength analyses. Nevertheless, 
the observed variability makes it very difficult to try to interpret any signal in 
terms of DM annihilation or decay in these sources. The Fermi satellite has 
more recently allowed the first detections in gamma-ray frequencies of nearby 
galactic objects, M31 [ITS], the SMC [175] and LMC [177]. As already em- 
phasized above, DM interpretations are difficult because of the presence of CR 
accelerators and gas in these objects, and these studies actually barely addressed 
the DM issue. 

Finally, many efforts have been recently done to check whether galaxy clus- 
ters could offer additional targets for DM searches (e.g. [1751 fTT^l \IM IT5T1 IT521 
183, 1184j ). Limits on DM candidates were again derived from the Fermi data 
[1851 11861 1181] , and we may also mention a recent claim for the detection of 
gamma-ray signals from a few clusters 187 . Nevertheless, wc know that high- 
energy processes are at work in galaxy clusters, and it is still difficult to infer 
and predict the related astrophysical background [1881 179J- Moreover, we stress 
that galaxy clusters host giant elliptical galaxies, like the M87 and NGC 1275 
galaxies discussed above. Since the latter have been demonstrated to be tran- 
sient TeV emitters, the relevance of galaxy clusters for DM searches might fairly 
be questioned, though their observations remain solidly motivated from many 
other astrophysical arguments. Deeper investigations appear to be necessary, 
which could for instance rely on the state-of-the-art cosmological simulations 
including baryon dynamics and star formation. 

7. Indirect dark matter searches with neutrinos 

High-energy neutrinos are interesting messengers to seek for DM annihilation 
traces. They can be produced through the fragmentation or decay of unstable 
primary annihilation products, or directly in some more rare cases. Direct pro- 
duction is generically suppressed if DM is made of Majorana fermions because 
of helicity suppression. Nevertheless, it may naturally appear in other cases, like 
in extra-dimensional models. Beside being complementary to other astrophys- 
ical messengers as tracers of DM annihilation in the Galactic halo, neutrinos 
do provide an additional and very specific signature: they can also trace the 
gravitational capture of DM in celestial bodies, in particular in the Sun or the 
Earth. Detecting high-energy neutrinos from the Sun or the Earth would rep- 
resent a very clean signature of DM capture and subsequent annihilation, hard 
to infer from other known phenomena. This idea was proposed in the mid 80's 
by several groups, notably by Press & Spergel |189j and Silk, Olive & Sred- 
nicki [190, for the Sun, and by Freese 191 , Krauss, Srednicki & Wilczek 192 , 
and Gaisser, Steigman & Tilav |193j for the Earth. 

Active neutrinos are almost massless and thereby travel in the skies like 
photons: their sources can be traced back, which is of particular interest for 
DM searches. Neutrinos are less subject to absorption than photons but in 
the meantime more difficult to collect on Earth. The detection of astrophysical 
neutrinos is indeed very challenging because of their weak cross section with 
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ordinary matter, which entails gigantic detectors, and because of the copious 
diffuse atmospheric neutrino background coming from interactions of cosmic 
rays with the atmosphere. So far, the history of neutrino astronomy relies on 
confirmed observations of two objects only: the Sun, from which the electron- 
to-muon neutrino flavor oscillation was proven, and the supernova SN-1987A, 
the explosion of which allowed to derive a robust upper limit on the absolute 
neutrino mass scale. These sources were observed in the keV-MeV energy range. 
At higher energies, more global (and indirect) astrophysical pieces of informa- 
tion have been delivered from e.g. measurements of the atmospheric neutrino 
anisotropy in the TeV energy range, related to the anisotropy in the CR flux. 

There arc currently two main types of so-called neutrino telescopes in opera- 
tion: MeV detectors historically optimized for solar/atmospheric electron/muon 
neutrinos — e.g. Super-Kamiokande^] (Super- K hereafter) — and GeV-PeV de- 
tectors mostly designed for muon neutrinos^] originating in high-energy astro- 
physical phenomena (supernova remnants, active galactic nuclei, DM annihila- 
tion, etc.) — e.g. AMANDA^IceCube^and ANTARES^KM3[^] These detec- 
tors are installed deep underground for protection against the CR background, 
and when dedicated to studies of astrophysical sources, the Earth itself is used to 
convert muon neutrinos into upward-going muons, the reconstructed directions 
of which allow a separation from the downward-going muons produced by inter- 
actions of cosmic rays in the atmosphere. Indeed, muons, when crossing such a 
detector, may produce Cherenkov light in the ice (AMANDA-IceCube are based 
at the South-Pole) or water (a tank of 50,000 tons of purified water for SuperK, 
the Mediterranean sea for ANTARES-KM3) that can be collected with arrays of 
photo- multipliers, allowing for direction reconstruction. The main background 
associated with this detection strategy comes from the atmospheric neutrinos 
generated by CR interactions in the atmosphere; although the amplitude of 
this background is large in the GeV-TeV range and makes neutrino astronomy 
challenging, predictions and measurements are available and under reasonable 
control (see e.g. |194|,ll95j ). Very approximately, the atmospheric muon neutrino 
flux is of order (j) hg (E) « 6x 10- 3 (£/100 GeV)" 37 m^s-isr- 1 . Assuming a full 
time exposure, an effective detection area of ~ 10 -2 m 2 (expected for a km-size 
neutrino telescope) and an optimistic fiat angular resolution of 1° [1961 1197] . 
this corresponds to 0(1 — 10) background event per year for a point-like source, 
among O(10 4 ) over half the sky. 

In the following, we introduce the treatment of neutrino propagation (in 
vacuum or/and matter) necessary for proper predictions in the present context. 
Then, we shortly discuss DM searches with high-energy neutrinos from the 
Galactic halo before delving into more details the case of solar neutrinos. We 



5. www-sk. icrr .u-tokyo . ac . jp/sk/ index- e .html 

6. Electron and tau neutrinos can also be reconstructed in some cases. 
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illustrate each paragraph with a selection of recent results and sketch the main 
prospects. 



7.1. Neutrino oscillations, interactions with matter, and muon flux at detectors 
Before discussing the Galactic or solar origin of the DM-induced neutrinos, 
it is interesting to review how they propagate to the Earth, and how one can 
predict the flux converted into muons at detectors. A fundamental feature of 
neutrino phenomenology is flavor oscillation, which is a consequence of their non- 
zero masses; interaction and mass eigenstates are different [198 . This of course 
must be accounted for in predictions, since high-energy neutrino telescopes are 
mostly designed for muon neutrinos. All cases are relevant here: oscillation 
in vacuum (Galactic halo neutrinos), to which one may further have to com- 
bine oscillation in matter (neutrinos from the Sun or the Earth). In the latter 
case, additional effects implying neutral and charge currents (absorption and 
re-emission) must also be plugged in. In the following, we ignore the potential 
effect of putative sterile neutrinos. We refer the reader to Refs. |199l 12001 1201) 
for extensive reviews. 

The neutrino oscillation framework (in vacuum) is usually based on a few 
ingredients: the mixing matrix U (including unknown but possible CP-violating 
phases), from which one can express a flavor (mass) eigenstate as a linear com- 
bination of the mass (flavor) eigenstates, and a classical quantum evolution 
equation assuming plane waves: 

\v a ) = U* ai \ Vi ), \v i {t))=e- iE * t \v i {Q))> (82) 

where Greek indices refer to the flavor basis (e, (i, t) while Latin indices refer 
to the mass basis (1,2,3). The oscillation probability between flavors a and /3 
is then merely given by V a p = \{vp\\i/ a )\ . Since masses are small, one can use 
the relativistic approximation such that + mf/2E, and 

the probability reads 

V aP (E, L = ct)= 6 a p-4 Tle{U* m U pi U aj U* pj ) sin 2 { ^ } (83) 

±2 £ lm(U* ai U 0i U aj U£ .) sin { ^ } . 

When distances get tremendous, e.g. when neutrinos originate in DM annihila- 
tion in the Galactic halo, the oscillation effects in vacuum average out and the 
probability is merely related to the raw mixing angles: 

(P a p(E)) L> ^f U 8 a p - 2 £ J>( Ue(U* ai Up U aj Ufa) . (84) 

Assuming no CP-violating phases, 023 — 7r/4 and #13 ~ 0, the neutrino flux at 
the Earth may be expressed as 202 

(85) 
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where any <f>° denotes the flux without oscillation. 

When neutrinos travel a significant distance inside a dense medium, the 
Mikheyev-Smirnov-Wolfenstein (MSW) effect [2031 1204] must also be included 
— this is relevant to neutrinos crossing the Sun or the Earth. Then, the above 
formalism has to be enriched with a new interacting term in the Hamiltonian 
which modifies the dispersion relation |2051 1206], and the flavor vector v — 
{u e , f M , v T Y now evolves according to 

v{t) = e- iHt u(0), (86) 

H = ^-Udiag^Aml^Am 2 ^)^ + 6iag(V2G F n e ,0,0), 
z hi 

where H is the Hamiltonian expressed in the flavor basis, Gf is the Fermi con- 
stant and n e is the electron density. The flavor asymmetry comes from that 
in coherent charged current interactions, electron neutrinos see the electrons 
present in the medium while the other neutrinos do not, leading to a flavor- 
dependent refraction index. One must also include other charged and neutral 
current interactions with nucleons (inelastic scattering), which get very impor- 
tant at high energy and mostly induce an energy redistribution toward lower 
energies because of successive absorption and re-emission processes (the con- 
version of tau neutrinos into tauons also induces extra flavor injection and a 
coupling in the evolution between neutrino and anti-neutrino species). This can 
be worked out with Monte Carlo methods |206j . An elegant and self-consistent 
alternative to such a two-fold treatment is that of the density matrix formal- 
ism |199( [207, 208J, which combines oscillation and absorption, and provides 
identical results. In that case, the flavor density matrix p v obeys the following 
evolution equation: 

dpu ,.dpu dp v dp v 

-tt = -t[H, p v + — + — + — , (87 

dl dl NC dl CC dl in 

where, in the right-hand side, the terms respectively describe the neutral cur- 
rent effects, the charged current effects, and the neutrino injection induced by 
DM annihilation. We emphasize that the injection term in a dense medium 
is different from what is obtained for annihilation in the Galactic halo, since 
some of the fragmentation or decay products of the primary particles, or even 
the primary particles themselves, can be absorbed before cascading. This is the 
case, for instance, for muons. Irrespective of the method, the full calculation 
has to be performed numerically. Furthermore, detailed models of the Sun or 
the Earth, when concerned, have to be used to get accurate predictions. 

An important outcome is that for neutrinos from the Sun, absorption is 
responsible for an exponential energy cut-off above 100 GeV in the differential 
flux. For neutrinos injected at the center of the Earth, the cut-off is somewhat 
higher, around 10 TeV. Moreover, it turns out that oscillation in the Sun mixes 
tau and muon neutrinos; oscillation only affects electron neutrinos in vacuum on 
their way to the Earth. This is is illustrated in Fig. [T3} which was extracted from 
Ref. |206j . In the top left panel, one can see the injected neutrino spectra at the 
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Sun's core, where the annihilation of 250 GeV WIMPs into t + t~ was assumed. 
In the top right panel, propagation effects (oscillation and weak currents effects) 
are shown by evolving the injected spectra from the Sun's core to its surface; 
the mixing between tauon and muon neutrinos is made clear, as well as energy 
redistribution toward low energies, due to weak scatterings; the production 
of secondary muon and electron neutrinos is due to tauon decays, which are 
generated from tau neutrinos through charged current interactions. Oscillation 
in vacuum between the Sun and the Earth (the distance is fixed to 1 AU) further 
results in slightly mixing electron neutrinos with the other flavors, as sketched 
in the bottom left panel. The final neutrino spectra at the detector are shown 
in the bottom right panel, where the oscillation pattern appears as completely 
smeared out. This does not come from additional mixing in the Earth, but 
from a time average performed over a 1 year cycle, during which the variation 
in the distance between the Earth and the Sun is enough to erase the oscillating 
features. Therefore, the net propagation effects for neutrinos coming from the 
Sun mostly translate into an amplitude rescaling, as in the case of Galactic 
neutrinos, and a high-energy damping due to efficient absorption in the Sun 
above 100 GeV. 

Once the neutrino flux at sea level is predicted, one still has to convert it 
in terms of muon flux at detectors (or other observables related to neutrinos, 



e.g. hadronic showers). Again, the full flavor evolution given e.g. in Eq. (87 1 
must be adapted to the Earth medium to determine the neutrino flux at a certain 
distance from the detector. For a muon-optimized detector, the relevant distance 
is the mean free path of a muon in the surrounding material, typically 1 km for 
a TeV muon travelling in the ice. Another important point to take into account 
is the time and detector-location dependence of the neutrino arrival direction. 
Indeed, only upward-going muons can be used for astrophysical searches, and 
detection efficiency is angular dependent. Therefore, the effective detection 
area A c s is a quantity somewhat tricky to extract for neutrino telescopes. It 
basically depends on the neutrino flavor, energy and direction, and has to be 
determined for each source. Anyway, for a given astrophysical source located 
at point x*, the detector efficiency might be time- averaged and the differential 
detected event rate then merely reads 



<av„(x») 

dt 



f dE„^- (A cS (E^)) . (88) 



The effective detection area is typically of order (*4 c ff) ~ 10~ 3 to 10~ 2 m 2 
(i?„/300 GeV) 2 for km-size instruments [196] . where one power of energy comes 
from the neutrino conversion into muon and the other from the muon energy 
loss converted into Cherenkov photons. To illustrate how challenging detect- 
ing astrophysical neutrinos is, this has to be compared with the usual effective 
areas encountered in gamma-ray astronomy. For instance, the Fermi-LAT in- 
strument^ has an effective area of ~ 0.5 — 1 m 2 in the 0.1-100 GeV energy 



11. http : //www-glast . stanf ord. edu/ 
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Figure 13: Evolution of the neutrino spectra between the Sun's core (top left), the Sun's surface 
(top right), the Earth's surface (bottom left), and the detector (bottom right). Legends are 
explicit, and more details are given in the text. These plots have been borrowed from Ref. |206l . 



range. 

7.2. Neutrinos from DM annihilation in the Galactic halo 

Since neutrinos are almost massless, they travel along geodesies like pho- 
tons, and may therefore trace the (squared) DM distribution in the Galaxy or 
extra-galactic objects. Thus, as for indirect DM searches with gamma-rays, the 
Galactic center on the one hand, and neighbor dwarf galaxies on the other hand, 
are very good targets. In these specific cases, the neutrino flux prediction can 
be written the same way as for gamma-rays, except for the possibility of flavor 
oscillations (see Sect. |7.lj ). For an arbitrary flavor index a, it reads: 
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<5Q bs t/los 



where V is the flavor oscillation probability, and p(l) is the DM mass density 
along the line-of-sight / in the chosen direction. The integral is performed over 
a cone of angle <5Sl bs corresponding to the angular resolution of the telescope 
(assumed flat for simplicity here) . This is therefore similar to the flux prediction 
made for gamma-rays except for the relevant angular size. Indeed, neutrino tele- 
scopes have generically angular resolutions worse than gamma-ray telescopes, 
ranging from a few degrees around a few GeV down to half a degree above a 
few TeV. Still, the predicted fluxes do not suffer from high-energy absorption 
like in the Sun, which is of interest for km-size neutrino telescopes, the effective 
detection areas of which typically increase like E 2 . 

Nevertheless, the predictions obtained for the Galactic center are generically 
small compared to the current experimental sensitivities, typically of the order 
of 0.1-1 up-going muon at detectors per km 2 and per year, much smaller than the 
annihilation signal expected from the Sun as it will be shown below (see e.g. 209, 
1210] ). Such predictions are usually based upon assuming a generic NFW DM 
profile, with an optimistic angular average of 1° around the Galactic center. 
The most optimistic estimate arises when requiring the potential complementary 
gamma-ray flux (model-dependent approach) to saturate the gamma-ray signal 
already measured at the Galactic center by different experiments, e.g. EGRET, 
HESS, Fermi. Proceeding so implicitly accounts for the still large uncertainties 
in our knowledge of the central DM distribution, which might be more cuspy 
than expected (e.g. because of the adiabatic compression triggered by baryon 
cooling), inducing fluxes significantly larger, around 100 up-going muons per 
km 2 and per year [209] . 

Conservative limits have been obtained for generic WIMPs fully annihilating 
into specific final states by the IceCube Collaboration |211) . They are reported 
in the (c ann v)-m x plane in Fig. [M] and shown to be rather far away from the 
relevant WIMP phase-space. This is essentially due to the poor exposure of 
IceCube to the Galactic center, since the experiment is located at the South Pole. 
Nevertheless, the analysis performed by the collaboration somehow compensates 
for this limitation by focusing on the high-latitude flux predictions, which are 
less subject to theoretical errors. Although this strongly damps the predictions, 
the DM density profile is in turn much less uncertain away than close to the 
Galactic center. 

Finally, we note that other interesting and similar DM sources are the neigh- 
bor dwarf galaxies, which are DM dominated objects. Nevertheless, expecta- 
tions are rather modest, and unless for TeV WIMP masses and decade timescale 
exposure, neutrino telescopes are poorly sensitive to such objects |213| . 

7.3. Neutrinos from DM annihilation in the Sun or the Earth 

The interaction properties of DM with ordinary matter open the possibil- 
ity of its gravitational capture and further accumulation in celestial bodies, in 
particular in the Earth and in the Sun which are the closest potential reser- 
voirs. Indeed, DM particles can lose kinetic energy through elastic scattering 
with ordinary matter and may thus remain gravitationally bound to the en- 
countered object — this happens when the final velocity is lower than the 
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Figure 14: Limits on the WIMP annihilation 
cross section obtained from observations of the 
Galactic halo at intermediate latitudes per- 
formed with the IceCube neutrino telescope. 
The unitarity bound derived in Ref. 121 21 is 
also shown. See the text for more details. This 
plot has been taken from Ref. [211] . 



so-called escape velocity. If sufficiently accumulated at the center, annihila- 
tion turns on and injects high-energy standard particles in the dense medium, 
among which neutrinos are the only species that can escape. As we will re- 
view below, this implies that the annihilation rate is not fixed by the (squared) 
average local DM density around the body, but is instead intimately linked to 
the capture rate. Once produced in the object, neutrinos experience the oscil- 
lation and absorption pattern discussed in the previous section before reaching 
the detector. We will focus on the Sun in the following, as a template ex- 
ample. For even more insights, we refer the reader to the seminal works in 
Refs. pEBHl [Ml QUI EEMJ CEH3 [2J1 [2T5], and also to subsequent developments 
found in e.g. Refs. [^1^1^1^1^122^123211205112051 1225112231 122511225] . 
Dedicated and sophisticated numerical codes can be found in Refs. |22ll 1206] . 

7.3.1. DM capture 

A pedagogical introduction to DM capture can be found in [227]. In the 
Sun, the total number N of DM particles evolves with time and depends on 
capture, annihilation, and escape, respectively as follows: 

d ^-=C c -C a N 2 -C e N, (90) 

although the latter can safely be neglected for particle masses above a few GeV 
|214]I215] . The capture rate T c is readily identified to parameter C c and depends 
on the elastic scattering cross section between DM and ordinary matter, while 
the annihilation rate reads 

r a (t) = V [ d 3 x(a ann «)n 2 (t,x) = °^N 2 , (91) 

J Sun 1 

where n is the (time-dependent) number density of DM particles. We will 
further explain how to evaluate these parameters. Before, let's have a closer 



look to Eq. (90). Interestingly, we note that when dynamical equilibrium is 



reached, i.e. dN x /dt — 0, the annihilation rate is then fully fixed by the capture 
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rate, with T a = C c /2. Nevertheless, this clearly depends on the annihilation 



parameter C a and the capture rate, and the general solution to Eq. ( 90 ) is given 
by: 

/ r< ( + ~\ . I r< 

(92) 




For neutrino detection today, this equation has to be evaluated at a time cor- 
responding to the age of the Sun, t — t Q = 4.6 Gyr. Dynamical equilibrium 
is reached when t@ j t c 3> 1, which is usually the case for most DM candidates. 
Indeed, typical values are of order t Ci © w 0.001 t c ^ ss 10 s yr — equilibrium is 
thereby barely reached in the Earth, while generically expected in the Sun. More 
precisely for the Sun, we have t Q /t c .q ~ y/ (C c /10 25 s _1 )((cr ann ^)/10 _30 cm 3 s _1 ), 
which means that even with an annihilation cross section much lower today than 
the one relevant at thermal decoupling in the early universe (this is the case when 
annihilation is p-wave dominated), the annihilation rate in the Sun can still be 



maximal. Indeed, in terms of annihilation rate, Eq. (92) gives 



r o (t) = ^tanh^j-j^^, (93) 

which is independent of the annihilation cross section when equilibrium is reached. 
We therefore anticipate that the solar neutrino flux can provide very strong con- 
straints to p-wave annihilation cross sections if the elastic scattering cross section 
is large enough, which is not the case with other indirect detection messengers. 
Given a point-like reservoir located at distance d from the Earth, the neutrino 
flux reads 

e v (E)T a dN v 
*» {E) = ^dT^E> (94) 

where dN v /dE is the neutrino spectrum injected at the source and e v (E) ac- 
counts for the oscillation and absorption effects discussed in the previous section, 
which can significantly alter the original spectrum. The point-like approxima- 
tion is usually very good for the Sun given the poor angular resolution achieved 
in neutrino astronomy, typically 1° or more for GeV-TeV neutrinos; this ceases 
to be the case for the Earth. In the following, we discuss the calculations of the 
annihilation parameter C a and the capture rate C c which both set the annihi- 
lation rate T a . 

To derive C a , one has first to guess the spatial distribution of DM in the 
Sun. Spherical symmetry applies, and if dynamical relaxation is assumed, then 
the spatial distribution is fixed by the local gravitational potential $ and the 
average kinetic energy (i.e. the average temperature T) |214) according to 

n(t,x) = n{t,r)=n {t)e- m x^ r '>/ T , (95) 

s f , , G N M(r') f r a 

$(r) = / dr n , M(r) = An I dr r p sun {r ) , 
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where no is the average density at the core, Cm is Newton's constant, M(r) 
is the solar mass within radius r, and p sun is the solar mass density. Since 
N — An J Q iisun dr r 2 n(t, r), it comes 



C a 



,u) n\t,r) 



N 2 



J^ m drr 2 e~ m x^)/T 



and we remark that we can express C a in terms of effective volumes as 



C a = 



V n = 



y'i 



2 ' 



(96) 



(97) 



3 TTlpj T^un 



-, 3/2 



2 n m x P; 



X Psun 



6.6 x 10^ 



lOGeV 



-3/2 



where the effective volume V n has an explicit dependence in the solar core 
temperature T sun and density p sun [227] . For the Earth, this volume is V n ~ 
2.3 x 10 25 (nrn x /10GeV) -3 / 2 cm 3 . Here, we have implicitly assumed that DM 
is concentrated at the Sun's core, and taken the corresponding constant values 
for the solar mass density and temperature. Instead, we may also evaluate 
these values at the mean orbital radius f [225J , so that we now have to solve the 
equation for C a using p sun (r) = p swa (f), an d T{r) = T(f). The mean orbital 
radius can be determined from the following implicit equation: 



d 3 x rn(t,r) _ dr r 3 e~ m * {f) 



-Rsun 



6T(f) 



J flsun d 3 x n(t, r) J Rsu ° dr r 2 e~ m * \ ^ 2 G N p sun (f) m> 



(98) 



This equation can be solved numerically. Then, the solution for C a is analytic 
and reads: 



C a 



V2 



(99) 



The authors of Ref. [225] have found that for masses larger than 1 GeV, a good 
numerical approximation to the previous result is given by 



C a ^ 4.5 x 10 



- 30 cm- 3 



0.6 GeV 



-i 3/2 



10 GeV 



(100) 



Here, the thermal average of the annihilation cross section (u ann v) should be 
performed according to a Maxwellian with the core temperature T(r) « 10 7 K 
fa 1 keV as the characteristic temperature. 

We now turn to the determination of the DM capture rate, which does 
depend on the scattering cross section and on the elemental distribution in the 
solar core |215j : 



c,0 



47T 
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Jo 



dr r 



\ - dTc,i( r ) 
^ dV 



(101) 
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It involves the differential capture rate per shell volume unit induced by a solar 
clement of label i: 



dT^r) _ Px Mr) *i f W - dw /H / x dQ FhQ) (102) 

Jo W J Q _ 



dV 2 m x p, 2 



where p x is the local DM mass density, p%{r) the element mass density in the 
spherical shell at radius r, at the WIMP-element scattering cross section, fa the 
WIMP-element reduced mass, f(w) the normalized velocity distribution in the 
solar frame, w the velocity as unaffected by the solar gravitational potential, 
F(Q) the nuclear form factor such that F(0) — 1, and Q the recoil energyp 2 ] 
The dependence on the elastic cross section and on the local WIMP density 
makes clear the correlation arising between the neutrino signal expected from 
the Sun/Earth and the event rate predicted at direct detection experiments, es- 
pecially when dynamical equilibrium is reached, in which case the annihilation 
rate is fully set by the capture rate. Additionally, we underline that given equiv- 
alent strengths for scalar (spin-independent) and axial-vector (spin-dependent) 
WIMP-nucleus interactions, capture in the Sun is usually dominated by the lat- 
ter because of the prominence of light elements in the core (mostly hydrogen), 
while the former plays an important role for capture in the Earth. 



In Eq. (102), the first integral is performed over the Galactic velocity w, 
which can easily be related to the velocity in the Sun v. Indeed, energy con- 
servation tells us that before the elastic scattering off a nucleus, the global 
WIMP kinetic energy oc m x w 2 remains constant. Besides, when the WIMP 
falls in the Sun's gravitational potential, gravitational energy gets locally con- 
verted into kinetic energy. Therefore, the WIMP velocity in the Sun is given 
by v 2 — w 2 + v 2 sc (r) , where the escape velocity v csc is defined as the velocity at 
which the kinetic energy equals the gravitational potential. The escape velocity 
at the Sun's surface is 617.8 km/s, to be compared to a value of 11.87 km/s for 
the Earth. A Maxwellian velocity distribution is usually assumed, which (in the 
Sun's frame) is given by 



f(w) _ / 3 1 



2 7T U v Vq 



3{w-v g ) 2 \ / 3(w + V 2 



exp < — -= > — exp 



2a 2 M 2ct 2 



(103) 



where Vq ~ 220 km/s is the velocity of the Sun in the Galactic frame, and 
a v ~ 270 km/s is the velocity dispersion of DM species. 

Kinematics must be inspected carefully, as detailed in |2151 1228. 23lj and 
very clearly summarized in 227 . A WIMP is captured as soon as its velocity 
is less than v esc after scattering. Defining 

(m x ± miY 



12. Form factors are different for scalar and spin-dependent interactions. Common refer- 
ences are [228112291 for the former and 1230) for the latter. 
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the fraction of WIMP energy lost in the collision, Jq = Q/E, must lie in the 
range 



< 

w 2 < 



2Q <(3+m x (w 2 +v 2 sc ) 



(105) 



for the WIMP to be captured. At this stage, it is interesting to remark that for 
an isotropic scattering (disregarding the form factor suppression for a moment), 
the probability V for a WIMP to scatter down to a velocity smaller than v esc 
should be flat in the above range such that [227 



V{w) 



1 



(w 2 



w 

T- 



w 



(106) 



where 9 is the Heaviside function. Incidentally, this equation shows that the 
probability of capture is maximized when /3_ is maximized, i.e. when m x rs 
rrii. This actually corresponds to resonance effects, which depend on both the 
elemental distribution in the celestial reservoir and the escape velocity. This 
has been intensely discussed in Ref. [228] where these resonance effects were 
predicted to be important for capture in the Earth, due to the larger relative 
density of oxygen and metals. 



Kinematics further helps define the ranges over which the integrals of Eq. ( 102 ) 
are performed: 
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(107) 



From these integral boundaries, we understand why a departure from the Maxwellian 
law can have a significant impact on the predictions, as discussed into more de- 
tails in e.g. Ref. |2"52] . 



7.3.2. Current limits on DM models and prospects 

We start by emphasizing that in contrast with other indirect detection mes- 
sengers, the solar neutrino signature is one of the best to discover DM annihila- 
tion because predictions do not suffer potentially big theoretical errors. Indeed, 
they mostly rely (i) on the knowledge of the local DM density and dynamics, 
which can be constrained reasonably well, (ii) on neutrino propagation in vac- 
uum and matter, which is under control, and (iii) on the DM particle mass, 
and its annihilation and interaction properties. The latter point can be fixed 
for each DM model. The most important uncertainty affects the evaluation of 
the WIMP-atom elastic scattering cross section, because of uncertainties in the 
estimates of the spin contents of nucleons (spin-dependent interactions) and of 
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Figure 15 : Top left : Muon flux predic- 
tions for supersymmetric DM models and 
limits from IceCube and Super-K. Bottom 
left: Corresponding limits in terms of spin- 
dependent WIMP-nucleon cross section. Bot- 
tom right: Corresponding limits in terms 
of spin-independent WIMP-nucleon cross sec- 
tion. 
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the nucleonic matrix elements (scalar interactions) — see e.g. [233] . Neverthe- 
less, there are hopes to reduce these uncertainties in the near future thanks to 
improvements in hadronic physics or QCD, for which experimental data will 
come from existing colliders. Finally, as another big advantage over other astro- 
physical messengers, the main neutrino background is well measured and also 
under control. 

Many predictions of the neutrino fluxes from the Sun and the Earth have 
been derived in the context of supersymmetric DM |227j . The reader can find 
extensive examples in e.g. |219l 1220, 2341 12351 1236] . Predictions in the frame 
of extra-dimensional theories are quantitatively similar (see e.g. |237l [238]). 
These predictions mostly concern the muon neutrino flux[^J which can be recon- 
structed at detectors as upward-going muons after charged current conversion. 
The predicted fluxes, usually higher if induced by DM annihilation in the Sun 
than in the Earth's interior, can reach 10 11 -10 12 upgoing muon neutrinos per 
year and square kilometer above a detection threshold energy of s» 10 GeV, 
which translates into an upgoing muon flux of 10 2 -10 3 km~ 2 yr _1 at detectors. 
The full energy dependence has further to be taken into account properly to 
convolve these expectations with the detector efficiencies and to compare them 
with the background. For DM searches toward the Sun, we note that in addition 
to the atmospheric background, a contribution is also expected from interactions 
of cosmic rays in the Sun itself, although it has been shown negligible |239|, 1240] . 

Disregarding the details of predictions, we emphasize that such a muon flux 
level is within reach of some current experiments like Super-K [241] and Ice- 



13. See Rcf. 226 for prospects with tau neutrinos. 
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Cube [242J , which have already derived interesting limits. These are shown in 
terms of muon flux in the top left panel of Fig. 15 taken from Ref. |242| . on 



top of predictions associated with DM candidates arising in the minimal su- 
persymmetric model (represented as the shaded area in the plot). Note that a 
broad range for the relic abundance was considered, down to a value as low as 
tt x h 2 = 0.05. This implies that the reported shaded area contains models with 
optimistically large cross sections, most of which corresponding to the largest 
predicted muon fluxes. Some direct detection constraints have been included, 
coming from both a combined analysis of the CDMS and EDELWEISS data 
|243| and from XENON 100 |244j . Two types of limits appear in the panel, one 
assuming a soft neutrino spectrum arising in annihilation into b quarks, and 
another one assuming a harder neutrino spectrum originating in annihilation 
into W bosons; the latter case is obviously more constraining. The supersym- 
metric models in the shaded zone, in contrast, usually come with more complex 
spectra, in between these two cases. While the derived limits are shown to start 
tickling the upper part of the relevant parameter space, the sensitivity of the 
complete version of IceCube is drawn as the red solid curve for a duty cycle of 
180 days. It clearly demonstrates the potential of neutrino telescopes to survey 
a large part of the supersymmetric parameter space. 

As already emphasized above, the neutrino signal induced by DM annihila- 
tion in celestial bodies should exhibit a correlation with direct detection signals. 
For the Sun, the neutrino signal is mostly proportional to the spin-dependent 
elastic scattering cross section, while the scalar one may also, in some case, 
become important. Therefore, limits on the solar neutrino flux can easily be 
translated in terms of elastic scattering cross sections and compared with the 
sensitivities of direct detection experiments — this is obviously DM model de- 
pendent. The bottom panels of Fig. [15] nicely illustrate that, where the muon 
flux predictions of the top panel are replotted in terms of WIMP-nucleon spin- 
dependent scattering cross section (bottom left) and of scalar scattering cross 
section (bottom right). We clearly see that, for the former case, current neutrino 
telescopes are already very competitive with respect to existing direct detection 
experiments like KIM [235] or COUPP [235]. On the other hand, the scalar 
cross section is more poorly constrained. 

The complementarity between direct DM detection and indirect detection is 
further illustrated in Fig. 16 taken from Ref. [225J , where the authors used the 
Super-K data to extract limits on generic WIMP models. The used data set 
slightly differs from the previous example in that only fully contained events, 
i.e. upgoing muons stopping within the detector, were considered to reduce the 
background and enhance sensitivity [224 . The top (bottom) panels show the 
limits obtained on the spin-dependent (scalar) cross section, while the left (right) 
panels concern models for which the relic abundance is set by a fully s-wave (p- 
wave) annihilation cross section. Note that a dominating p-wave annihilation 
cross section makes indirect DM detection with other messengers irrelevant since 
annihilation in galactic halos is then velocity suppressed; neutrinos are there- 
fore the only indirect detection possibility in this case. In the plots, the authors 
have assumed different pure annihilation final states, which are made explicit. 
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Figure 16: Limits on the spin-dependent (top) and spin-independent (bottom) WIMP-nucleon 
cross sections obtained from the Super-K data. The left (right) panels assume that the WIMP 
relic density is set by a pure s-wave (p-wave) annihilation cross section. Regions favored by 
some direct detection experiments appear as shaded areas. These plots have been taken from 
Ref. [225] , 



Finally, they have reported as shaded zones the regions corresponding to the 
excess events found in two direct detection experiments, DAMA |247j and Co- 
GeNT [37], together with the limits obtained by the PICASSO [MB], COUPP 
[MS], and CDMS [MS] collaborations. 
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